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TRAJECTORY MODULE OF THE NASA AMES RESEARCH CENTER 

Aircraft synthesis program acsynt 

Michael E, Tauber and John A. Paterson - ' 
Ames Research Center 


SUMMARY 


A program was developed to calculate trajectories for both military and 
commercial aircraft for use in the NASA Ames Research Center aircraft synthe- 
sis program, ACSYNT. The function of the trajectory module was to calculate 
the changes in the vehicle's flight conditions and weight, as fuel is consumed, 
during the flying of one or more missions. The traj ectory calculations 
started with a takeoff, followed by up to 12 phases chosen from among the fol- 
lowing: climb, cruise, acceleration, combat, loiter, descent, and paths. In 

addition, a balanced field length was computed. The emphasis was on rela- 
tively simple formulations and analytic expressions suitable for rapid compu- 
tation since a prescribed trajectory had to be calculated many times in the 
process of converging an aircraft design, or finding an optimum configuration. 
The trajectory module consists of about 2500 cards and was operational, as 
part of the ACSYNT program, on a CDC 7600 computer. 


INTRODUCTION 


Computerized aircraft design at the conceptual, or early preliminary 
design level, has proven to be a practical engineering tool (refs. 1-3). The 
development of such a program requires the efficient integration of many dis- 
ciplines such as the aircraft's geometric layout, weight estimation, aerody- 
namics, propulsion, and trajectory analysis. This report contains a descrip- 
tion of a program, which was developed to calculate trajectories for both 
military and commercial aircraft and was used in the NASA Ames Research Center 
aircraft synthesis program, ACSYNT. The function of the trajectory module was 
to calculate the changes in the vehicle's flight conditions and weight, as 
fuel is consumed, during the flying of one or more missions. The trajectory 
calculations started with a takeoff, followed by up to 12 phases chosen from 
among the following: climb, cruise, acceleration, combat, loiter, descent, 

and paths. In addition, a balanced field length was computed. The emphasis 
was on relatively simple formulations and analytic expressions suitable for 
rapid computation since a prescribed trajectory had to be calculated many 
times in the process of converging an aircraft design, or arriving at an opti- 
mum configuration. The trajectory module consists of about 2500 cards and was 
operational, as part of the ACSYNT program, on a CDC 7600 computer. Although 
it could potentially be modified, the program is presently not Suitable for . 
stand-alone computations since it depends on inputs from a propulsion module 
for thrust aiid fuel flow rate, and an aerodynamics module for lift and drag 
coefficients. 
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GENERAL PROGRAM! DESCRIPTION 


The function of the trajectory module is to calculate the changes in the 
flight conditions and weight that take place in a vehicle during the flying of 
one or more missions. Any number of missions may be flown during a given com- 
puter run, where each mission is divided into phases, and particular phases 
may be further subdivided into legs. A given mission may have up to 12 phases, 
selected in any order from among the following 7 types: climb, cruise, accel- 

eration, combat, loiter, descent, and paths. The trajectory module (called 
TRAJ hereafter) consists of 16 different subroutines organized as shown in the 
block diagram below. The arrows point in the direction of called subroutines. 


FROM ACSYNT CONTROL 



A brief description of the above subroutines follows; 

TRAJA - This is the master control program for TRAJ. ICALC is a control 

variable passed from ACSYNT control to TRAJA and informs TRAJA what 
tasks are to be performed. 

TRAJIN - Reads all input data for TRAJ. 

TRAJOI - Prints all input data for TRAJ. 

MISSON — Secondary control program for TRAJ. Calls all necessary subroutines 
needed to fly one mission. Subroutines for each phase of the mission 
are called in the same order as the various phases of the mission are 
read in by TRAJIN. 
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TAKEOF - Performs all calculations required during takeoff. This is always 

the first routine called by MISSON. It. is not counted as a phase of 
the mission. 

CLIMB — Performs all calculations of a climb phase. 

CRUISE - Performs all calculations of a cruise phase. 

ACCEL - Performs all calculations of an acceleration phase. 

COMBAT — Performs all calculations of a combat phase. 

LOITER - Performs all calculations of a loiter phase. 


DESCNT - Performs all calculations of a descent phase. 

PATHS — Performs all calculations of a paths phase. Input consists of up to 
50 pairs of altitude vs velocity (or Mach number) data points. 


! .... 

GOLDEN - Performs golden section type of search for (1) Breguet altitude or 
! (2) maximum endurance altitude or (3) optimum endurance Mach number. 

f 

: CONVGR - Converges vehicle on gross weight. Process consists of taking esti- 

mated gross weight, calling on MISSON to get computed gross weight, 
comparing estimated and computed gross weight to get a new value of 
estimated gross weight, and repeating this process until vehicle 
weight has converged to within a specified tolerance limit. 


TRAJ00 - Prints all output data for one mission. 

TRPLOT - Computes and writes onto tape various parameters which are to be 
displayed by some graphics device. 


CONTROL PROGRAM - TRAJA 


This is the master control program for TRAJ. ACSTNT control communicates 
with TRAJ only through this program. When ACSYNT control calls TRAJA, it 
always sends a value for a variable called ICALC which tells TRAJA what tasks 
are to be performed by TRAJ. 

Briefly, when ICALC - 1, TRAJ reads and prints all TRAJ input data. When 
ICALC = 2, TRAJ executes for mission 1 (the basic mission) . When ICALC = 3, 

TRAJ makes one final execution of mission 1 and prints output. Further, if 
other missions exist (missions 2, 3, 4, . . .), then these missions are exe- 
cuted oyer aiid over until they have converged, on. gross "takeoff weight,- with the 

airframe kept fixed as it existed after mission JL was completed. The output 
of these converged missions is then printed.'-- * • 

After mission 1 has been completed, the fixed airframe weight is computed 
by the equation i 
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WIXED = WGTO - WTOT - TOL 

where WGTO is the gross takeoff weight of the converged vehicle, WETOT is 
the total fuel weight, and WL is the payload weight. 

INPUT-TRAJIN 

This routine reads all the TRAJ input data. The data consist of two 
parts. The first part consists of general input data which is not specific to 
any one mission. This part is stored in core memory only. The second part is 
the data relevant to some particular mission. This second part is stored on 
scratch tape where it can be read back into core memory, mission by mission, 
as it is needed. 

The general 7RA.J input data consists of the following 40 variables: 

(1) TIMT01 - Time during takeoff at first power setting 

(2) TIMT02 - Time during takeoff at second power setting 

(3) FREURE - Fuel reserve, as fraction of total fuel 

(4) DESLF - Design load factor 

(5) ULTLF - Ultimate load factor 

(6) RANGE - Aircraft range 

(7) WFUEL - Total internal fuel weight 

(8) WFEXT - External fuel weight 
. (9) JJETRAP - Weight of internally trapped fuel 

(10) MENDUR - Mach number used for final endurance calculation 

(11) QMAX - Maximum dynamic pressure 

(12) XDESC - Initial guess for horizontal distance covered during descent 

(13) WKFUEL - Fuel multiplication factor 

(14) CRMACH - Cruise Mach number 

(15) NCRUSE - Number of cruise legs in mission 

(16) IFSIZE - Power setting for sizing engine 

(17) IPST01 - Power setting during first part of takeoff 
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(18) IFST02 - Power setting during second part of takeoff 

(19) ' IBREG - Breguet altitude search indicator 

(20) IENDUR - Endurance altitude search indicator 

(21) - IPRINT - Diagnostic print indicator 

(22) KERR.0R. - Error print indicator 

(23) NLEGCL - Number of legs per climb phase 

(24) MLEGCR - Number of legs per cruise phase 

(25) NLE6L0 - Number of legs per loiter phase 

(26) MILCOM - Takeoff obstacle height indicator 
The following quantities are used in the landing field computation: 

(27) WKLAND - WLAND = WGTO - m<\HD * WTOT 

(28) ETtFAC - FLLAND = (XAIR + XGRLAN ) /FLFAC 

(29) DEGEL - Deceleration factor 

The following quantities are used for the plotting routine: 

(30) IPLQT - Plotting indicator 

(31) HMINP - Minimum altitude 

(32) HMAXP - Maximum altitude 

(33) DELHP - Altitude interval to be used 

(34) SMMINP - Minimum Mach number 

(35) SMMAXP - Maximum Mach number 

(36) DEIMP - Mach number interval to be used 

(37) WCOMBP - Aircraft weight used 

The following quantities are used for multimission computations: 

(38) NMISS - Number of missions per job 

(39) FWGMAX - Maximum Weight factor used in convergence 

(40) TOL - Convergence tolerance 


The input for each mission consists of 2 variables which apply to the 
whole -missions and a series of 14 variables which appxy to each phase of the 
mission. 

The two general variables are: 

(!) NPHA.SE - The number of phases in the mission 

(2) /PL - Payload weight (not needed fpr mission 1, where it is supplied 
by the weights module) 

The following 14 variables are input for each phase of the mission, 
although no single phase makes use of all 14 variables . 

(1) MSTART - Starting Mach number. If MSTART > 0, then starting Mach number 

is set equal to MSTART. If MSTART < 0, then starting Mach 
number is set equal to final Mach number from preceding phase. 

(2) MEND - Ending Mach number 

(3) HSTART — Starting altitude. If HSTART > 0, then starting altitude is set 

equal to HSTART. If HSTART < 0, the starting altitude is set 
equal to final altitude from preceding phase. 

(4) HEND — Ending altitude. In CLIMB, if IBREG “ 0 or HEND > 0, then 

ending altitude is set equal to HEND; otherwise, ending altitude 
is determined by making a : Breguet search. 

(5) X - Horizontal distance covered 

(6) TIM -Time 

(7) NT - Number of turns 

(8) VIND - Constant indicated airspeed (for climb) 

(9) IP - Power setting 

=1 maximum afterburner power 

-.2 intermediate (commercial maximum takeoff power) 

=3 maximum continuous power 
=4 thrust = drag (etuis e) 

» 5 idle 

(10) IX - Indicator which tells whether horizontal distance covered during 

CLIMB, ACCEL, or DESCNT will be added, subtracted, or ignored 
in computing horizontal distance to be covered during CRUISE 
phase.'. 

(11) It? - Weapons drag indicator 

(12) IPRT - Print indicator 
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(13) IB - Bombs drop indicator - 

(14) IA AMMO drop indicator 

The table below indicates which inputs may possibly be needed for each of 
the seven types of phases. 


. PHASE. 

1 

MSTART 

2 

MEND 

3 

HSTART 

4 

HEND 

5 

X 

6 

TIM 

1 

NT 

8 

VIND 

9 

IP 

10 

IX 

11 

1W 

12 

IPRT 

13 

IB 

14 

IA 

CLIMB 

/ 

v 

V 

/ 




/ 

V 

/ 

/ 

/ 

/ 

✓ 

CRUISE 

/ 


/ 


/ 




/ 

/ 

/ 

/ 

/ 

/ 

ACCEL 

/ 

/ , 

/ 
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/ 

/ 

/ 

/ 

/ 

COMBAT ' 

V 


/ 



/ 

/ 


/ 


/ 

/ 

/ 

/ 

■ ■.LOITER 

: : / 


/ 


/ 

/ 





/ 

/ 

/ 

/ 

DESCNT 




/ 







/ 

/ 

✓ 

/ 

PATHS . 

Mach 

vs h table 






/ 


/ 

/ 

/ 
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MISSION 


The primary function of this subroutine is to call all the subroutines 
needed to fly one mission. For example, if a given mission consists of three 
phases such as climb, cruise, and descent, then MISSON will first call TAKEOF 
(always called first for all missions), then CLIMB, then CRUISE, then DESCNT. 
The order of the called subroutines is determined by the order of the various 
phases of the mission as read by TRAJIN. 

After this, MISSON computes used, reserve, total, and internal fuel 
weights according to the following formulas: 

WFUSED = WKFUEL * (WGTO - W - f(W)) 

WFRES = WFUSED/ (1 . - FRFURE) - WFUSED 
WFTOT = WFUSED 4- WFRES + WFTRAP 
WFUEL « WFTOT - WFEXT 

where WFUSED is the weight of the fuel used, WGTO is gross takeoff weight, 

W is the landing weight of the vehicle, f(W) is the weight of all missiles, 
bombs, ammunition, and external tanks jettisoned during the mission, WKFUEL is 
a multiplying factor, WFRES is the reserve fuel weight, FRFURE is the fraction 
of fuel kept in reserve, WFTOT is the total fuel weight, WFTRAP is the weight 
of fuel trapped in the fuel tanks, WFUEL is the weight of internal fuel, and 
WFEXT is the weight of fuel stared in external tanks. 

Next, the landing field length is determined by the following sequence of 
equations (ref . 4) i 

WLAND - WGTO - WKLAND*WFUSED 
VSTALL = SQRT [ 2 . *WLAND / ( SWING* CL*RHO ) ) 

VSCREN = 1* 3* VSTALL 


VTDOWN - 1.15*VSTALL 

XA2R - [ (VSCREN**2 - VTDOWN)/ 64- 4 + 50.]*RLD 
XGRLAN = VTDOWN**2/(64.4*DECEL) 

FLLAND = (XAIR + XGRLAN) /FLFAC 

where WKLAND, DECEL, and FLFAC are input to TRAJ (see TRAJIN writeup), RHO is 
sea level air density, CL is lift coefficient at landing, and RLD is lift to 
drag ratio. 

Next, the endurance altitude and endurance time are calculated based on 
the approximate (input) Mach number (MENDUR). (The results of the approximate 
calculation for endurance altitude and time were only intended for use in the 
economics module of ACSYNT and are never used within TRAJ for performance or 
fuel calculations.) A golden section type of search is performed to find the 
endurance altitude corresponding to the maximum value of the endurance factor 
(RLD/SFC). Using this value of the endurance factor, the endurance time is 
calculated from RLD/ SFC [WGTO/ (WTO - WFUSED) ] . 

Finally, the gross weight of the vehicle is calculated using the equation 

WGCALC - WFIXED + WPL + WFTOT 

where WPL is the payload and WFIXED is equal to WGTO - WFTOT - WPL calculated 
just once after mission 1 has been converged. WGCALC is needed for transfer 
back to CONVGR. 


Takeoff 


The first calculation performed in subroutine TAKEOF is the determination 
of the total static thrust of the aircraft, and the corresponding value of 
SFC. The fuel used during startup and taxi is based on the static value of the 
SFC for a power setting which is input (XFST01) and for an input time 
(TIMT01). 


The takeoff balanced field length is computed using an expression from 
reference 4, which is written 


FLTO = 2.10 


0 ,01163 W 

PCt S 

‘-"7 


+ 0.374 TOOBHT 


] 


+ 2.7 + 


32.0 


°‘ 04 


where FLTO is in ft, W/S in psf and p in slugs/ft 3 . The lift coefficient 
C L is evaluated at 1.2 V gtal ^, where 


stall 




2 WGTO 

pC T S 
L Stall 


k t is the thrust lapse factor during takeoff and the program uses the average 
of the value at takeoff and the value at a velocity of 1.2 "00BHT 
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is the obstacle height and is 50 ft (15.24 m) for military aircraft and 35 ft 
(10.67 m) for commercial aircraft. 

The fuel used during takeoff can be calculated by one of two methods . If 
the takeoff time (TIMT02) is input as a positive or zero value, the fuel used 
iS 

TJFT02 = 0. 0167 (TIMT02) (SEC) (TN) 


where the SFC and thrust (TN) are based on conditions at the end of takeoff. 
Alternatively, if TIMT02 is set negative, the fuel used is based on the bal- 
anced field length and is 


■RFT02 


«W <SF W 

1800. (1.2 V a£aU ) 


where and are values averaged between static and takeoff 

conditions. 


The total fuel used is the sum of the values for startup , taxi and 
takeoff 


WTO = WT01 + WTO 2 


Climb 

The program can calculate climb in either of two modes. When VIND > 0, 
climb is calculated using a constant indicated airspeed. When VIND < 0, 
climb is calculated using an approximate minimum fuel path. If the aircraft 
is to cruise at Breguet altitude (IBREG > 0}* the program performs a search 
for the ending altitude; otherwise ending altitude must be input (IBREG = 0). 

The forces acting on the airplane during climb vary with altitude, neces- 
sitating dividing the climb phase into a finite number of legs. The condi- 
tions changing with altitude include the temperature, pressure, and density of 
the ambient air s and the weight and Mach number of the vehicle. Also, since 
the available power decreases at an ever greater rate as the ceiling altitude 
of the vehicle is approached, it is desirable to decrease the altitude incre- 
ments of the legs as altitude increases. 

The total number of legs in a climb phase is given by the formula 


NLEGS - -EgtgL friterval 

average step size 













f where Hj is the starting altitude of the phase, is the ending altitude 

of the phase, and f(H) is an expression giving step size as a function of 
altitude. A further restriction on NLEGS is given by the relation 
3 < NLEGS < 20. The above formula can be overridden by giving NLEGCL a posi- 
tive integer /alue, in which case, NLEGS = NLEGCL. This, however, will cause 
all climb phases in all missions of a run to have the same value for NLEGS. 

Once the number of legs in a climb phase is determined, the distance 
climbed during each leg is given by the formula: 


■l 





DELH = 


1 + (LEG) (e - 1) / N LEGS ] 

1 + (LEG - l)(e - 1)/NLEGS_P H 2 


Hj) 


where LEG is the current leg being executed by the program* This causes the 
step size to decrease logarithmically with increasing altitude* 


Summing forces in the axial direction, the equations of motion give 

/T D \ 

aF = - w - sin V 

while in the normal direction 


L + T sin ;* = W cos ) 

ttfhich gives two equations in three unknowns. The three unknowns are T, L , an< 
D, or v * 


For the approximation to the most economical or fastest climb, the rela- 
tion for y comes from reference 5, and is 


3 

I e - + V (s)V + 1 

where 
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For the constant indicated airspeed climb. 




where, now 


V. , V. , 
xnd _ rnd 

^P/PSL ^ 


and where 0 is the density ratio referenced to the sea level value* By dif- 
ferentiating the equation for velocity with respect to time, we find that the 
relationship between time and altitude is 


At 


/a .. V ind 
Ah - ■ 


V. , 
md 


2go 


1.5 


T p SL V ind C D S 
W " 2W 


We now have all the necessary terms for calculating flight-path angle. 
The fuel used is calculated from 

EW_, = Z (sf c)T At 

r 


and the equivalent ground distance travelled during climb is 

EX - EV cos y At 


Cruise 

The starting altitude for the cruise phase will be the input value 
(HSTART) if HSTART is positive* If HSTART = 0 and IBREG « 1, then the 
starting altitude will be determined by a Breguet search* -The ending altitude 
mil be the same as the starting altitude only if the input value (HEftl)) is 
set equal to -1* Otherwise, the airplane will be allowed to climb with each 
leg of the cruise phase as fuel is burned off. 

At tne beginning and at the end of the cruise phase a test is made to 
determine whether the fuel used to that point exceeds the weight of the fuel 
stored in the external tanks if these are present. If it does, then the 
external fuel tanks are dropped immediately. 

In general, the horizontal distance travelled in a given cruise phase is 
calculated using the formula 

XCRUSE = X + BXCRUS + XDESC 

where XCRUSE is the actual horizontal distance that will be traversed in the 
cruise phase, X is the nominal horizontal distance input for that phase, 
DXCRUS is the algebraic sum of the horizontal distances traversed in the pre- 
ceding climb and acceleration phases (but following the preceding cruise 
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phase if one exists), and XDESC is the horizontal distance traversed in the 
last descent phase of the mission. An estimated value for XDESC is inputted 
and used for the first iteration of the mission, and the calculated value from 
the previous iteration is used for all succeeding iterations of the mission. 
DXCRUS is set to zero at the beginning of the mission and is reset to zero at 
the end of every cruise phase. The calculated value of horizontal distance 
traversed in a climb or acceleration phase will be added or subtracted from 
DXCRUS or ignored depending on whether IX for that phase is +1, -1, or 0, 
respectively. Similarly, XDESC will be positive, negative, or zero depending 
upon whether XX for the cruise phase under consideration is +1, -1, or 0, 
respectively. 

When X is given a negative value, then the formula for XCRUSE becomes 
XCRUSE = RANGE/NCRUSE + DXCRUS + XDESC 
where RANGE and NCRUSE are input at the beginning of the run. 

The number of legs in a given cruise phase (NLEGS) is given by the fol- 
lowing inequalities: 

X < 100, then NLEGS « 2 
100 < X < 300, then NLEGS - 3 
300 < X < 1000, then NLEGS = 4 
1000 < X < 3000, then NLEGS = 5 
X > 3000, then NLEGS = 6 

The cruise consists of an approximation to the mpst efficient flight path 
which uses a stepwise cruise-climb path. There are two options for determining 
the initial cruising altitude: (1) the Breguet altitude, found by maximizing 

the value of (V/sfc)(L/D) as i function of altitude and (2) altitude is 
specified. 

The summation of forces tangential to and normal to the flight path are, 
respectively. 


T cos a = D and L = "W 
From these relations, the weight of fuel used is 




C W(sfc) f — — - 

cos a W Ax 


where AX is the distance flown at each altitude and the cruise altitude is 
increased, as the aircraft burns fuel, according to the relation 
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and C L is kept, constant at- the value corresponding to (L/D) ina j : ; otherwise, 
is determined by the input values of Mach number and altitude. 

The cruising time and the equivalent ground distance covered are, 
respectively, 

t = Z — and X - EAX 


Acceleration 

The purpose of the acceleration phase (ACCEL) is to increase the Mach 
number of the vehicle from an initial low value to a specified higher value. 
The phase is divided into steps with the size of the Mach number steps being 
0.02 between Mach 0.87 and 1.1 and equal to 0.05 elsewhere. 

The aircraft is assumed to accelerate at a fixed altitude. The equation 
of motion along the flight path is 


dV /T D\ 

dt 8 \W C0S a “ W/ 


from which we can calculate the fuel used to be 


W F - h|i - exp[§J" if 


(sfc)T 
cos a - D 


The time required to accelerate is 




WdV 

T cos a - D 


and the equivalent ground distance covered is given by 


-iJ 


WdV 

T cos ct - D 


Combat 

In the combat phase, the acceleration and turning capabilities are calcu- 
lated an<’ also the amount of fuel used by performing combat at a fixed alti- 
tude, ej *'*r for a specified time or a fixed number of turns. The airplane’s 
potential auceleiation capability is expressed as the "specific excess power," 
for which the expression is derived below. 
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The specific energy (sum of potential and kinetic energy per unit air- 
plane weight) is 

E V 2 

E S = I " h + 2^ 

and differentiating the specific energy with respect to time gives the spe- 
cific power 

p = fk _ dh V dV 
s dt dt g dt 

Summing forces along the flight path gives 

dV 

T cos a ~ D - W sin y = m 

Dividing both sides by W and multiplying by V 

V(T cos a - D) „ . , V dV 

— — = V sm y + — -rr 

W g dt 

and, since 


dh 

dt 


- V sin y 


we get 


dF 

p _ s _ V(T cos a - D) 

s “ dt W 


The relations for turning radius and turning rate are, respectively 
(ref. 6) 

V 2 d0 • V 

R = — r r and -rr = 0 = pr 

g tan <J> dt R 

where the airplane bank angle, cj>, is related to the load factor, n, through 

, 1 

COS <p = — 

n 

The load factor is written as 


n = 


L 4- T sm a 
W 


and cannot exceed the structural design load factor. 
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Two types of turns are considered. These are: (1) sustained, or coordi- 

nated turns for which P s = 0, squiring the solution of the equation 
T cos a « D to obtain the aerodynamic parameters needed to find the load 
factor and (2) instantaneous turns corresponding to the maximum value of P s 
and limited by either the maximum lift coefficient, or the structural 

design load factor of the airplane. 

The time required to complete a turn is 

_ 2ttR 
~ V 

and the fuel used during combat turns is 

W F = ET(sfc)Nt 

where N is the number of turns required. 


Loiter 

The starting altitude for the loiter phase will be the input value 
(HSTART) if HSTART is positive. If HSTART = 0 and IENDUR = 1, then a 
search will be made to determine the most economical starting altitude. 

The Mach number will change from leg to leg or remain constant throughout 
the loiter phase depending on the input value of MSliRT. Ordinarily, a search 
for the most economical loiter Mach number will be made for each leg of the 
phase. However, whenever MSTART is positive, the Mach number will remain 
equal to MSTART throughout the phase and the altitude in each leg will 
increase to take advantage of the reduced fuel weight. 

Each loiter phase is divided into legs, the number of which depends on 
the amount of time to be spent in the phase. The number of legs is equal to 
the time in hours rounded up to the next higher integer, but is not allowed 
to exceed six. 

There are two options for determining the loiter flight conditions: 

(1) most economical loiter altitude and/or velocity is found by maximizing the 
value of (L/D)/sfc as a function of altitude, velocity, or both, and (b) alti- 
tude and/or velocity is specified. 

The fuel used is calculated from 

XW_ = ET (sf c) At 
r 

where t is the loiter time. As the aircraft burns fuel, the loiter altitude 
is increased, or the speed decreased, to maintain the most economical flight 
condition according to the relation 
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Since loitering generally consists of flying in circles, it may be 
desirable to fly in a smaller circle than that corresponding to the most 
economical one. To decrease the loiter circle size, more lift is needed. The 
equation of motion normal to the flight path is 


L + D tan a - 


1 + & 


where R is the loiter circle radius. The exact solution is a tedious itera- 
tion between C^, Cq, and a, since 


^ " C D tan a = ^ 


1+1 — 

\Rg J 


However, for most cases of practical interest 


C D tan a << C L 


and, we have used the simple approximation 



Descent 

The most economical descent path is the flattest descent with engines at 
idle power setting. The aircraft should fly at (L/D) m a v , thus = Cp o and 

for flight with thrust much less than drag, the optimum descent angle is 
approximately (ref. 6) 


sin y = - 


Summation of forces normal to the flight path gives 


L = W cos v 


resulting in the following expression for the speed 

.. /2W cos y 


The descent time is calculated from 


EAt « E 


V sin y 


while the fuel used is 

EW_ = ET (sfc)At 
I 1 

and the equivalent ground distance covered is 

LAX - EV cos Y At 

OUTPUT- TRAJOO 

The output quantities are listed in tabular form and consist of 19 quan- 
tities representing values at the end of each phase. In addition, there are 
eight values giving a fuel weight. If there are combat phases, an additional 
21 values of combat parameters are printed for each combat phase. Lastly, 
there are 12 miscellaneous quantities printed. 

The following 19 quantities are printed for each phase: 

(1) Mach number 

(2) Altitude 

(3) Lift coefficient 

(4) Angle of attack 

(5) Fuel used in phase 

(6) Flight time in phase 

(7) Flight velocity 

(8) Specific fuel consumption (installed) 

(9) Specific fuel consumption (uninstalled) 

(10) Total thrust (installed) 

(11) Total thrust (uninstalled) 

(12) Drag coefficient 

(13) Lift-drag ratio 

(14) Weight at end of phase 

(15) Engine airflow rate 

(16) Dynamic pressure 

(17) Horizontal distance covered 



(18) Engine installation drag coefficient 

(19) Pressure recovery 

The following eight quantities give a fuel weight: 

(1) Fuel weight during takeoff x*ith first power setting 

(2) Fue*l weight during takeoff with second power setting 

(3) Mission fuel weight 

(4) Reserve fuel weight 

(5) Trapped fuel weight 

(6) Internal fuel weight 

(7) External fuel weight 

(8) Total fuel weight 

The following seven quantities are printed for each combat phase for each 
of three conditions: lg flight, sustained, and instantaneous — making a total 
of 21 combat parameters in all. They are: 

(1) Specific excess power 

(2) Load factor 

(3) Turning rate 

(4) Radius 

(5) Angle of attack 

(6) Lift coefficient 

(7) Drag coefficient 

The following 12 miscellaneous quantities are also printed out: 

(1) Total mission time 

(2) Total mission range 

(3) Takeoff field length (total run) 

(4) Landing field length (total run) 

(5) Landing field length (ground run) 





(6) Weight used for landing calculation 

(7) Takeoff weight 

(8) Landing weight 

(9) Endurance Mach number (input value) 

(10) Endurance altitude 

(11) Endurance time 

(12) Loiter radius 


PLOTS-TRPLOT 

This routine calculates various parameters which are written onto tape 
so that they can later be displayed graphically. The main part of this rou- 
tine calculates certain parameters which may be displayed graphically on a 
plot of altitude vs Mach number. These parameters are; 

(1) PSIG = specific excess power at Ig 

(2) NZS = load factor during turn 

(3) TDOTS = sustained turning rate 

(4) NZI = maximum instantaneous load factor 

(5) PSI = maximum instantaneous specific excess power 

(6) TDOTI = maximum instantaneous turning rate 

(7) Q = dynamic pressure 

These parameters are calculated for a certain range of values of altitude and 
Mach number with beginning, ending, and incremental values of each being 
input . 

The inputs to the routine are; 

(1) IPLOT = 0, no print, no plot 

= 1, print, no plot 
= 2, no print, plot 
= 3, print, plot 

(2) HMINP — minimum altitude 
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(3) UMAX? = maximum altitude 

(4) DELHP = altitude interval 

(5) SMMINP = minimum Mach number 

(6) SMMAXP = maximum Mach number 

(7) DELMP = Mach number interval 

(8) WCOMBP £ 0, aircraft weight taken from first COMBAT phase 

0 < WCOMBP < 1, aircraft weight equals WCOMBP times takeoff weight 
WCOMBP > 1, aircraft weight equals WCOMBP 

The second part of this routine calculates landing field lengths at alti- 
tudes of 0, 4000 ft (1219 m), 8000 ft (2438 m) , and 12,000 ft (3657 m) . 



PROGRAM LISTING 
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Zt 


25 


3 C 


35 


AC 


A 5 


5 C 


55 


TRAJA 76/75 OPT -2 


FTN A ■ 5 + A 1 Q l 


SUBROUTINE TkAJA l IC ALC , NERROR/ LGEO/ KG PRNT , I G PLT 1 

EXTERNAL MIS5QN 

REAL MSTART/NEND/N7/MENDUR 

CCMMON /TRAJCM/ ALPHA/ ARM* CD/CDL/ CDO/ CL/ DE SL F / DR AG# EN# HN/ KP/ 

1 R LD/SKN/QMAX /RANGE/ SFCz SWING/ THRUST /TWzULTLF/W/WAMMUN, WE TANK/ 

2 WFU£L/WGT0WT/WMISS/WT0T/WPLWT/MSTART(1Z)/HSTART(12)/TIM( 12)/ 

3 PS1GT(12)/TD0TST(12)/NZST(12)/PSIT (12 ) / TDOT I T C 12 ) / NZI T { 1 2 ) / 

A WBOMBS/ WF EXT/ SFCU/THRSTU/ COINS P/ PR TOT / WKFUEL/ C RMACH/ F LTO / X( 12) j 
5 FLL AND/ TENDURz SLRANG/ BLTIMEz TIMTOl /TIMTQ2/WFT0/PNAMEKX2)/ 
t PNAMEZt L2)/SMNTC12)/HNT(12)/WFm2)/TINtT(12)/XTtl2 ),CLTC12)/ 

7 CDT(12)/ALPHATC12)/CLIT(1Z)/CDIT(12)/ALPHIT(12)/WFUSED/WLAND/ 

8 um2>/RLDT(l?)/$FCT(12)/TNT(12)/CET( 12)/0Y (18)/ 

9 IAO/IPS/ITS/IHS/IBS/NPHASE/IDT (1A) 

COMMON /TRAJEX/ AL P1GT <12 ) /COINS T(12 ) / CD1GT 1 12 J / C LlGT C 12 ) / DEC EL/ 

1 DELHP»DELMP/0XCRUS/FLFAC/FRFURE/HENDtl2)/HENDUR/HMAXP/HMINP/ 

2 HPATHT (53)/MEND( 12 ) / HENDUR/ NT { 12 > / PRTOTT £ 12 ) / RAD 1 T C 12 ) / 

3 RADST{12)#SF/SFCUT(12)»SMHAXP/SMMINPz 

A TNUT(12)zVIN0<12)/VELTCl2)/VPATHT(50)/WART£12)zWC0MBP/WEELT(12)z 
5 WFRES/WFTOT/ WF TQ1/ WFTQ2/ WFTRAP / WKL AND / XDESC/ XGRL AN/ F WGMA X/T OL/ 
b WFIXED/WGCALC/WPL/WGT0/DUMY(82)/IA( 12)/ IAS/ I B ( 12 ) / I BREG/ IENDUR/ 

7 IMISS/ IPU2J/IPHASE/IPL0T/IPRINT/ IPRT (12) /IP SIZE / IPSTQl/ IPST02/ 

8 IW(12)/IX(12)/KERROR/MILCOM/NCRUSE,NLEGCL/NLEGCR/NLEGLO/ NMISS/ 

9 NPA TH/ I DUKY ( 6 ) 

DATA PATH/AHPATH/ 

WPL-WPLWT 
WGTO=*WG TOWT 

IF (KGPRNT»EQ.Q) GO TO 5 

KSTORE-KERROR 

KERR OR *2 


5 

IF 

{ IC ALC • c C • 1 ) 

GG 

70 10 


IF 

( ICALC.EG. 2 ) 

GO 

TO 20 


IF 

( 1 CALC.E 0 . 3 ) 

GQ 

TO 30 


IF 

( KGPRN T * Nt . 0 

1 KERROR-KSTORE 


F^LTURN 



C— 


— 



c 

ICALC 

- 1 



c — 







10 TIMT 01 * 5.0 

TIMTU 2 - 1.0 
FRFURE-.C 5 
DESLF- 2.5 
LLTLF- 3.75 
RANGE- 1000 . 
InFUE L* 10 jOC . 
FFLX T* 0«0 
WF 7 RAP* 100 . 
CMAX* 700 . 
XUESC- 80 . 
WKFUEL* 1*0 
LkMACH-.B 
WHLAND -.57 
FLFAC -.6 
DECEL «.25 
NCRUSE *2 
IPS IZE *0 
1 PST 01*5 
IPST 02*2 


ORIGINAL PAGE IS 
OP POOR QUALITY 
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SUBROUTINE TRAJA 


76/76 OPT-2 


FTN 4.5+4 10 


C 


6G 


65 


70 


75 


85 


9C 


95 


10C 


105 



IBREG-0 

I ENDUR *0 

IPRINT-0 

KEKROR-2 

I PLOT* J 

HMINPO.O 

HHAXP-40000. 

DELHP-4000. 

SMKINP* .3 

SHHAXP-.9 

DELHP-.l 

WC0MBP-.7 

MLCOM-O 

MENDUR* .5 

NLEGCL-0 

NLEGCK-C 

NLEGLD-0 

NhI5S»i 

FkGMAX* 1 . 2 

TOL=.C0l 

CALL TRAJIN 

CALL TRAJGI 

C READ HISS I ON 1 FROM TAPE 2Z. 

12 IMISS-1 
KEhtlND 22 

READ ( 22/ 22C0 ) NPHASE/ WPLTRA 

2200 FCJRMaT(I10/E15.8) 

DC IS IPHASE-l/NPHASE 

READ ( 22/2201 ) PNAKE1 ( XPHASE >/ PNAME2 ft PHASE ), MSTARH IPHAS E )/ 

1 MENU(IPHASE)/HSTART(IPHASE)/HEND(IPHASE)/X( I PHASE }> TIHCI PHASE 

2 NT ( IPH ASeJ/VINDt IPHAS E) / IP ( IPHASE ) / IX ( IPHAS E )/ IW( IPHASE) / 

3 I PR T( IPHASE) / ISdPHASE ) / I A t IPH AS E } 

2201 FORM AT ( 2 A4/2X/7E15. 3/1:15. 8/6110 ) 

IF (PNAMEl(IPHASE) .NE.PATH) GO TO 15 

READ (22/2200) NPATH 

READ (22/2202) ( HP ATHT (I) / 1-1/ NPATH ) 

READ (22/2202) ( VPATHT ( I ) / 1-1/ NP ATH ) 

2202 FORMAT ( 8E15. 8 ) 

15 CONTINUE 

C — ---END OF READING HISSIQN 1 FROM TAPE 22. 

IF (ICALC.EC.3) RETURN 

WF IXED«0. 0 

SF-ULTLF/OESLF 

ALPHA«0.0 

CL *0.0 

DRAG-0.0 

HN-0.0 

SNN-0.0 

h -0. 0 

IF (KGPRNT.NE.O) KERRQR-KSTORE 
RtTURN 

C 

C ICALC ■ 2 

20 CALL HI SSONC ICALC/ NERRQR/ IGEO# KGPRNT) 

IF (KGPRNT.NE.O) KERRGR-KSTORE 
RETURN 
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SUBROUTINE TRAJA 


76/76 


OPT" 2 


FTN 4.5+4, 


0 


115 


120 


125 


130 


135 


14C 


145 


150 


155 


16C 




c 

C ICALC - 3 

30 IF (NHISS.EQ.l) GOTO 40 
KCALC-2 

CALL MlSSQN(KCALC,NERROR,IGEO,KGPRNT) 

IF (NERR0R.GE.2) RETURN 

4C CALL TRAJOU 

IF (IHISS.GE.NMISS) GO TO 70 
IF (IMISS.GE.Z) GO TO 50 
WGE5T-WGT0 
WGMAX-FWGMAX+WGEST 
VFIXED"WGTO— WFTOT-WPL 

C READ NEW MISSION FROM TAPE 22. 

50 IMISS-IMISS+1 

READ (22,2200) NPHASE,WPLTRA 
DO 60 IPHASE*1,NPHASE 

READ (22,2201) PNAM61 ( IPH ASE ) >PNAME 2 t I PHAS E ) , HST ART ( IPH AS E ) , 

1 MEND ( I PHASE ), HST ART ( I PHASE ) , HENO ( I PHA SE ) , X t IPHASE ), TIM ( 1 PHASE), 

2 NT(IPHASE),V.'ND(IPHASE),IP(IPHASE>, IX { IPHAS E ) , I W ( IP HASE) , 

3 IB ( IPHASE 1 , I A ( IPHASE) 

IF (PNAMEUIPHASE).NE.PATH) GO TO 60 

READ (22,2200) NPATH 

READ ( 22, 22C2 ) t HP ATHT ( I ) » I»l# NPATH ) 

READ (22,2202) (.VP ATHT ( I ),.I.;J,,.NPATH.) 

60 CONTINUE 

C END OF READING NEW MISSION FROM TAPE 22. 

WPL-WPLTRA 

WGEST1-WGEST 

IF (IPRINT.EG.O) GO TO 65 

WRITE (6,600) I MISS, WGEST, WGHAX , WF IXEO ,WFTGT,WGTO#WGCONV,WPL, 

1 WGEST 1, WGCALC 

600 FORMAT ( 1H1, 1 BH CHECK FOR MISSION, I2/23H BEFORE CALLING CONVGR*/ 

1 9H WGEST*, E13.6,4X,7H WGMAX", E13 • 6, AX, 7HWF IXEO", E13 . 6/ 

2 9H WFT0T»,E13.6,4X,7H WGTQ«,E13.6> 4X, 7HWGC0NV", E 13.6/ 

3 9H WPL",E13.6,4X,7HWGEST1«,E13.6,4X,7HWGCALC«,E13.6) 

65 CALL CONVGR(MISSON,WGEST1,WGCONV, WGHAX, TOL,IPRINT) 

IF (IPRINT.EQ.O) GO TO 40 

WRITE (6,601) WGEST, WGMAX, WF I XED, WF TOT, WGTO, WGCONV, WP L, WGE ST1, 

1 WGCALC 

601 FORMAT ( H2H AFTER CALLING CONVGRi/ 

1 9H WGEST*, E13.6,4X,7H WGHAX", E13. 6, 4X, 7HWF I XED", E 13 . 6/ 

2 9H WFT0T",E13.6,4X,7H WGTO", E13 . 6, 4X, 7HWGC0NV-, E 13 . 6/ 

3 9H WPL",E13 t 6,4X,7HWGESTl",E13.6,4X,7HwGCALC«,E13.6) 

GO TO 40 

70 IF ( I PLOT ,GT .. 0) CALL TRPLDT( NERROR, IGEO,KGPRNT) 

IF (KGPRNT oNC.O) KERROR"KSTORE 

GOTO 12 

END 
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10 


15 


20 


25 


30 


36 


40 


46 


50 
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SUBROUTINE TRAJIN 

REAL MSTART/ MEND* NT/ MENOUR 

CGMMQN /TRAJCM/ ALPHA/ ARW/ CD/ CO L/ CDD/ CL# DE SL F » DR AG* ENz HN# KP/ 

1 RLD/SMN/UMAX/RANGE/SFC/SWING/ THRUS T/ T W/ UL TL F / Wf WAWMUN/ WE TANK* 

2 WFUEl/WGT0WTzWMISSzWT0T/WPH/T/M5TART(12)/H5TART(12)/TIMI12)/ 

3 PSlGTtl2)/T0Qrsr(J.2IfN£5T(l2)zPSIT(12)fT00TIT(lZ)zNZIT(12J/ 

4 W8DMBS/ WFEXT / SFCU/ THR3TU#CDINSPfPRTQTfKKFUELfCRMACH/FlT0fXll2l/ 

5 FLLANO/TENDUR/8LRANG/BLTIME/TIMTG1/TIKTOZ'WFTG/PNAME1(12)/ 

6 PNAKE 2 U 2 )zSMNTU 2 )/HNT<lZ)zWr'T{lZ)f riMET‘i.2)/ XT< 12 )/CLTU 21 » 

7 CDT(12)/ALPHAT(12 > / CL IT ( 12 ) / CD I T( 12 ) / ALPNIT 112 )/ WFUSED/WLAND/ 

8 QTtl2)zRL0T(12)»SFCT(12l/TNTtl2J/CET(12)/DY(18}/ 

9 IA0/IPS/ITS/IWS/IBS/NPHA3E/ IDY(14) 

COMMON /TRAJEX/ ALP1GT l 12 ) z CDINST (12 ) / CD1GT (12 ) / CL1GT ( 12 J / DECE It 

1 OELHP/OELMP/ DXCRUS/ F LFAC/ FRFUREz HENO ( 12 i / HE NOUR* HMAXP/HM INPz 

2 HPATHT(50)/MEND(12)/MENDUR/NTtl2)zPRT0TT(12)>RADITtl2)# 

3 RA0ST(12)/SF»SFCUT(12)/SMMAXP/SMHINPz 

4 TNUT(12)/VlND(12)/VELT(lZ)/VPATHT( 50) / WARTC12 )/ WCOMBP / HE 6 >, T 1 12 ) / 

5 WFRES/WFTUT/WFTOl/WFTGZzWFTRAP/WKLANO/XDESC/XGRLAN/FWGMAX/TQl/ 

6 WFI XED/ WGC ALC/ WPL/ WGTO/OUMY ( 82 ) # I A { 12 ) / I AS/ IB ( 12 )z IBREG# IENDURf 

7 I MISS/ IP l 12) /IPHASE/ IPLOT/ 1 PRINT/ I PRT 1 12 ) / 1 PS I ZE / 1 P STDlf I PSTOEj. 

6 IW(12J»IX(12 l/KERRCIft/MILCaMzNCRUSEzNUGCLzNLbGCRzNLEGLOz NMISSz 
9 NPATH/ IDUMYtfc) 

DATA PATH/4HP ATM/ 

NAHEUST/TRDATA/riMT01/TIHT02/FRFUREf DESLF/ULTLFfRANGE/WFUELf 

1 WFEXT/WFTRAP/QMAX/XDESC/ WKFUEL/ CRMACH/WK.LAND/FLFAC/DECEL./NCRUSE/ 

2 IPSIZE/IPSTOlz IPST02/ IBREGt IENDURz I PR INT/KERRORz IPLOT/HMINP/ 

3 HHAXPf DELHP/ SHHINPfSHHAXPfDELHPfWCOHBPf HIuTOlif HENDURf NLEGCLy 

4 NLEGCR/NLEGLQ/NHISSzFHGMAX/TGL 
READ 15/TRDATA) 

REWIND It 

00 10 IhISS b 1/MHISS 

READ (5/500 MPHASE/ WPLTRA 

WRITE ( 22» 2200) NPHASE/ WPLTRA 

500 FORMA T( I10/E10.0) 

2200 FQRMAT(I1Q/£15.8) 

C READ IN THE FOLLOWING WORDS FOR EACH PHASE 

C P H A S F. ( 2 PARTS)/ STARTING HACH/ ENDING MACH/ STARTING ALTITUDE/ 

C ENDING ALTITUDE/ HORIZONTAL DISTANCE/ TIME/ NUMBER OF TURNS/ 

C CONSTANT INDICATED AIRSPEED CLIMB/ POWER SETTING/ RANGE INCREMENT 
C INDICATOR/ WEAPONS RELEASE INDICATOR/ BOMBS RELEASE INDICATOR/ 

C AMMO RELEASE INDICATOR. 

DO 10 IPHASt*l/NPHASE 

READ (5/501) PNAME1 { IPHASE ) / PNAHE 2 1 1 PHAS E ) » MS TAR T ( IP HASE > / 

1 MENDt IPHASE) /HSTARTt IPHASE) /HENO (IPHASE )/X( IPHASE) /TIMl I PHASE)/ 

2 NT ( IPHASE)/ VINOC IPHASE) /IP (IPHASE) #IX{ IPHASE) jlWt IPHASE)/ 

3 IPRTt IPHASE )>IB(IPHASE)/IA( IPHASE) 

501 FORMAT ( 2A4/2X/7E 10.0/ E 10. 0/6110) 

WRITE (22/P2G1) PN AME 1( I PHAS E ) / PNAME2 ( I PHASE ) / MS T ART (I PHA S E 1 * 

1 MENDt IPHASE )/HSTART( IPHASE) /HEND( I PHASE )/X( IPHASE) /TIM (I PHASE)/ 

2 NT( IPHASE )/ VINO ( IPHASE )/IP( IPHASE)/ IX ( IPHASE )/IW( IPHASE) » 

3 IPRTt I PHA St )/ I B( IPHASE )/IAt IPHASE) 

22G1 F0RMAT(2A4,2X/7E15.8/E15.6/6I10J 

IF (PNAMEl(IPHASE).NE.PATH) GO TO 10 

READ (5/500) NPATH 

READ (5/502) ( HP ATHT ( I ) / I«lz NPATH ) 

READ (5/502) ( VPATHT ( I ) / I «1/ NPATH ) 

5G2 FORMAT ( 8EL0 .0 ) 
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WRITE t 22* 2200 ) 

WRITE l 22* 2202 ) 

WRITE (22*2202) 

f ORMAT ( 8E1 5 • 8 ) 

CONTINUE 

RETURN 

END 


NPATH 

(HPATHT(I)*I-1»NPATH) 

(VPATHTII)>I*1*NPATH) 
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1 SUBROUTINE TRAJOI • 

C PRINTOUT OF TRAJECTORY INPUT 
REAL MSTART,MENO,NT,MENDUR 

COMMON /TRAJCM/ ALPHA, ARW,CD,CDL,CDO,CL,D ESL F> DR AG, EN,HN,KP, 

5 1 RLD,SMN,QMAX, RANGE, SFC, SWING, THRUST, TW, ULTL F, W> W AMMUN, WE T ANK, 

2 WFUEL, WGTQWT ,WKISS,WT0T,WPLWT,MSTART(I2},HSTART(12),TIM(12), 

3 P51GT(12),TDOTST(12),NZST(12),P$1T(12),TDOTIT(12),NZIT(12), 

4 WBDMBS,WFEXT,SFCU,THRSTU,CDIHSP,PRT0r,WKFU£L,CRMACH,FLT0,XU2), 

5 FLLAND,TENDUR, BLR ANG, BLTIHE, TIMT01,TIMTQ2,WFT0, PNAMEH12), 

It 6 PNAME2(12),SMNTE12),HNT(12),WFT(12),TIMET(12),XT(12),CLT(1Z), 

7 CDT(12),ALPHAT(121,CLIT{12>,CDIT112), ALPHIT ( 12 ) , WFUSED, WL AND, 

8 UT(12),RL0T(12),SFCT(12),TNT(12),CET(12),0Y(18), 

9 I AO, IP S, ITS, IW$,IBS,NPHASE,IDY(14) 

COMMON /TRAJEX/ AL P1GT( 12 ) , CDINST( 12 ) , CD1GT < 12 ) , CL1GT( 12) , DECEL, 

15 1 DELHP,DELMP,DXCRUS, FLF4C,FRFURE,HEND(12),HEN0UR,HHAXP»HMINP, 

2 HPATHT(50),MEND(12), ME NDUR, NT( 12) , PRTOTT 1 12 ) , R AD IT ( 12 ) , 

3 RADST(12),SF,SFCUT(12),SMMAXP,SMMINP, 

4 TNUT(12),V1ND(12),VELT(12),VPATHT(50),WART(12),WC0HBP,WEELT(12), 

5 WFRES,WFTOT,HFI01,WFXP2iWFTRAP,i/KLAND,XOESC,XGRLAN,FWGMAX,TOL, 

2C fc V«FIXED,WGCALC,WPL,WGTO, DUMY(B2),IAC12),IAS, IB(12 ), I8REG, IENDUR, 

7 IMISS,IP( 12) , IPHASE, IPLOT, IPRINT, I PRT { 12) , I PSIZE, IPST01, I PST02, 

B IW(12),IX(12) , K ERROR , MIL COM, NC RUSE, NL EGCL, NLEGC R, NL EGL 0, N MISS, 

9 NP ATH, IDUMY (6 ) 

DATA PATH/4HPATH/ 

25 WRITE (6,600) TIHT01, KENDUR, NCRUSE, I PL OT, TIMT02, QMAX, I PSI ZE, HMINP, 

1 FRFURE,XDESC, I PST01, HMAXP, DESL F, WKFUE L, I PST02, DELHP , ULTL F, 

2 CRMACH,IBRE G, SMMINP, RANGE, WKLAND, I ENDUR, SMM AXP, WFUE L, FLF AC, 

3 IPRINT,DELMP,WFEXT,DECEL,KERROR, WCOMBP, WFTR AP, NL EGC L, NLE GCR, 

4 NL£GLO,FWGMAX,TOL»MILCOM,NMISS 

30 600 F0RMAT(1H1,25X,16HTRAJECT0RY INPUT/// 

1 IX, 8HTIMT01 «,F4.1, 7X, 8HMENDUR -,F6.0, 4X,BHNCRUSE *,I2, 

2 3X,8H IPLOT *, 12 / 1X,8HTIMT02 -,F4.1, 7X,8H QMAX *,F7.0, 

3 3X, 8H IPS IZE *, 12 , 3X,BH HMINP *,F7.G/ IX, 8HFRFUR E «>F5.Z, 

4 6X,8H XDESC *,F6.1, 4X,8HIPST01 «,12 , 3X/8H HHAXP «,F7.0/ 

35 5 IX, 8H DESLF »,F5.2, 6X,8HWKFUEL -,F6.3, 4X,8HIPST02 *,I2, 

6 3X,8H DELHP *,F7.0/ IX, 0H ULTLF »,F5.2, 6X,8HCRNACH *,F6.3, 

7 4X,8H IBREG «,I2 , 3X, BHSNMINP ",F6.3/ IX, BH RANGE *,F7.0, 

8 4X, 8HWKL AND *,F6.3# 4X,8HIENDUR -,I2 , 3X,8HSHMAXP «,F6.3/ 

9 IX, 0H WFUEL «,F8.0, 3X,8H FLFAC *,F6.3, 4X,8HIPRINT «,I2, 

4C 1 3X,8H DtLMP «>F6.3/ IX, 8H WFEXT «,F8.0, 3X,8H DECEL «,F6.3» 

2 4X» 0HKERRQR «,I2 , 3X,8HWCUHBP *,F9.2/ 1X»0HWFTRAP -,FS.O, 

3 3X, 8HNLEGCL >,I2 , 8X, 0HNLEGC R -,I2 , 3X, 0HNLEGL 0 *,12/ 

4 IX, 8HFWGMAX «,F6.3, 5X,8H TQL *»F6.3, 4X,8HHILC0M *,I2, 

5 3X,8H NKISS «,I2) 

45 REWIND 22 

DC 10 IMISS«1,NMISS 
READ (22,2200) NPHASt,WPLTRA 
2200 F0RMAT(I10,E15.8> 

IF (IMISS.EQ.l) WRITE 16,601) IHISS 
50 601 FORMAT ( / / / 24X,7HMISSI0N,I2// 

1 9X, 66HMACH MACH ALT ALT HORIZ NO. IP 

2 IX IW/75H PHASE START END START END DIST TIME TURNS 

3 VIND IPRT IB IA/) 

IF CIMISS.GT.1) WRITE (6,602) IMISS,WPLTRA 
55 602 FORMAT ( /// 24X, 7HMISS ION, 1 2, 12H (PAYLOAD *,F8.0,4H LB)// 

1 9X, 66HMACH MACH ALT ALT HORIZ NO. IP 

2 IX IW/75H PHASE START END START END DIST TIME TURNS 
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3 VINO IPRT 10 IA/) 

DO 1C IPHASE*!, NPHASE 

READ (22,2201) PN AMEK I PH ASE ) , PN AHE2 ( I PHASE > , MST ART ( I PHAS E ) , 

1 MfcNO(IPHASE),HSTARTl IPHASE > , HEND UPH A Sc ) , X ( IPHASE I, TIM (I PHASE), 

2 NT (IPHASE 5 , VINO (IPHASE I, IP (IPHASE ), XX (IPHASE ),IW( I PHASE) , 

3 IPRTdPHA St), IB (IPHASE), I At IPHASE) 

2201 FORM AT (2A4,2X,7E15.8/£15.0, 6110) 

WRITE ( 6, 6C3 ) PNaMEKIPHA5E),PNAME2(IPHASE),MSTART(IPHASE) , 

1 HENO ( I PH ASEKHSI ART (IPHASE), HE ND(I PHASE ),X(1PHASE),TIM( IPHASE), 

2 NT( IPHASE »> VINO (I PHASE KIP (IPHASE), I X( I PKASE),IW( I PHASE) > 

3 IPRTt IPHASE), IB ( I PHASE), I A( IPHASE) 

603 FORMAT (IX, A4, A3, F 5. 2, F 6 . 2, 2F8 .0 , F 7. 1, F 6. 1, F6 . 1,F7. 1, 14, 15, 14/ 

1 165,15,14) 

IF (PNAKEKIPHASEJ.NE.PATH) GO TO 10 

READ ( 22, 22C0) NPATH 

READ (22,2202) (HPATHT(I), 1*1, NPATH) 

READ (22,2202) (VPATHT(I), 1*1, NPATH) 

22C2 FORMAT (BE 15. 8) 

WRITE (6,604) NPATH 

WRITE (6,605) (HPATHT(I), 1*1, NPATH) 

WRITE (6,606) (VPATHTII), 1*1, NPATH) 

WRITE (6,607) 

604 FORMAT ( 9X» 8H NPATH -»I3> 

6C5 FORMAT ( 9X, 6HHP ATHT *, 7Fd . 0/ ( 17X, 7F8 .0) ) 

6C6 F0RMAT(9X,8HVPATHT *, 7F8 . 0/ ( 17X, 7F8 . 0) ) 

607 FORMAT! /) 

10 CONTINUE 
RETURN 
END 
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SUBROUTINE hIS SON 76/76 QPT-2 FTN 4.5+410 05 


1 SUBROUTINE MISSONI IC ALC# NERROR# IGEO#KGPRNT) 

REAL LQIT»MENDUR 

COMMON /TRAJCM/ ALPHA*ARW#CD#CDL#CDO#CL#DESLF#DRAG#EN#HN#KP# 

1 R LO#SHN# QHAX# RANGE# SFC #S WING# THRUS7#TW#ULTLF# W# WAHHUN# WE TANK# 

5 2 WFUEL# WGTGWT# WMLS S# W TOT# WPLWT# MS TART (12)#HSTART(12) #TIH(12) # 

3 PS1GT(12)#TDQTST(12)#NZST(12)#PSIT(12)# TDOTIT <12 )#NZIT (121# 

4 wBOMBS#WFEXT#SFCU#7HRSTU#CDINSP#PRTOr#WKFUEL#CRMACH#FLTQ#X(12)# 

5 FLLAND tTENDUR# BLR ANG#BLTIMZ#TI HTQ1# TI HT02# WFTQ#PNAME1 C 12 ) # 

6 PNAH£2{123#SHNT(12i#HNT(12)#wFTU2>#TIHET<12>»xni21#CLT(lZ)# 

10 7 CDT(12)#ALPHAT(12)#CLIT(12)#CDIT(12I#ALPHIT(12)# WFUSEO#WLAND# 

B QT(12)#RL0m2).SFCT(12t#TNT(l2)>CET(12 5#DY(18)# 

9 IAO#IPS#ITS»lWS#IBS#NPHASE#IDY 114) 

COMMON /TRAJEX/ ALP1GT 1 12 ) # COINS TU2 ) # CD1GTI12) #CL1GT C 12) # DECE L# 

1 DELHP#DELMP#DXCRU5#FLFAC#FRFURE#HEND(12)#HENDUR#HMAXP#HMINP# 

15 2 hPATHT(50)#MEN0(12)#MENDUR#NT(12)#PRTaTT{12>#RA0IT(12)# 

3 RADST(12)#SF#SFCUT(12>#SMMAXP#SMMINP# 

4 TNUT( 12 )#VIN 0 ( 12 )#l/ELT( 12 )#VPATHn 50 ) # WARTU2 >» WCOMBP# WEELT 112) # 

5 WFRES»WFT0T#WFT01#WFT02#WFTRAP#WKLAND#X0ESC#XGRLAN#FWGHAX#T0L# 

6 WFIXED#WGCALC# WPL# WGTQ# DUMY ( B2 )# I A 1 12 J # I AS# IB 1 12 ) # I BRSG# IENDUR# 

2C 7 IMISS»IPtl2)#iPHASE#IPLQT#IPRINT#IPRT(12 ) #IPSIZE# IPSTQ1# IP5T02# 

b lKtl2}#IX(12)#KERRQR#MILC0M#NCRUSE#NLEGCL#NLEGCR#NLEGL0#NMISS# 

9 NPATH# I0UMYI6) 

DATA CLIM/4HCLIM/ 

DATA CRUI/4HCRUI/ 

25 DATA ACCE/4HACCE/ 

DATA C0MB/4HCDMB/ 

DATA L0IT/4HL0IT/ 

DATA DESC/4HDESC/ 

DATA PATH/ 4RPATH/ 

3C C SET ULTLF SO AS TO MAINTAIN RATIO BETWEEN ULTLF AND DESLF 

20 ULTLF-S F+DESLF 

C SET WEAPONS# BOMBS# AMMO# AND TANKS DROP INDICATORS TO 1. 

C THIS WILL CAUSE AERO TO CALCULATE DRAG INCREMENTS 
C FOR WEAPONS# BOMBS# AND TANKS UNTIL SOME TRAJECTORY PHASE 
35 C SETS THEM TU ZERO. IF THE AMMO INDICATOR IS CHANGEO 

C TO ZERO BY SOME TRAJECTORY PHASE# THIS WILL AFFECT 
C THE CALCULATION OF WFUSED. 

IWS-1 

IBS-1 

4C IAS-1 

IT5-1 

C SET SUH OF HORIZONTAL DISTANCES COVERED DURING CLIMB AND 
C ACCELERATION PHASES TO ZERO AT BEGINNING OF MISSION. NEW VALUE 
C WILL BE CALCULATED IN’EACH CLIMB AND ACCELERATION PHASE AND 
45 C WILL BE SET BACK TO ZERO AT END OF EACH CRUISE PHASE. 

DXCRUS-0,0 

C SET BLOCK TIME AMD BLOCK RANGE TO ZERO SO THAT THEY CAN BE 
C ACCUMULATED IN EACH PHASE. 

BLTIME-0.0 
5C BLRANG-O. 0 

CALL TAKEOF ( ICALC# NERROR# IGEO# KGPRNT) 

IF CNERRQR. GE. 2) GO TO 160 
IF (W.LT.WPL) GO TO 13Q 

C PERFORM ALL PHASES OF MISSION IN THE ORDER THEY WERE READ IN. 

55 DD 120 IPHASE-1#NPHASE 

IF ( PNAME1 ( I PHASE ) .EQ.CL IM) GO TO 40 
IF (PNAHEl(IPHASE).eQ.CRyi) GO TO 50 
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lit 


601 

40 

50 

60 

70 

80 

90 

100 

11C 

120 
C — 
130 


602 

C — 
14C 


C 


rF (PNAHEKIPHASE) .EQ.ACCE ) GQ TO 60 
IF (PNAHEKIPHASE). EO. COMB) GO TO 70 
IF (PNAHEKIPHASE). EG. LOIT) GO TO 80 
IF (PNAHEKIPHASE) .EQ.DESC) GQ TQ 90 
IF (PNAHEKIPHASE), EQ. PATH) GO TO 100 
WRITE ( 6* 601 ) IPHASE 

FORHAT (/47H TRAJECTORY INPUT USES UNDEFINED 
GG TO 160 


CALL CLIKB( IC ALC* NERRQR* IGEQ* KGPRNT ) 
GU TQ 110 


NAME FOR PHASE* 13/ ) 


CALL CRUISE ( IC ALC* NERRQR* IGEO* KGPRNT) 

GO TO 110 

CALL ACCELt ICALC*NERROR* IGEO.*KQPRNT} 

GQ TQ 110 

CALL COMBAT (ICALC/NERRQR* IGEO*KGPRNT) 

GO TQ 110 

CALL LOITER ( ICALC*NERROR* IGEO^KGPRNT) 

GQ TQ 110 

CALL OESCNKICALC* NER RQR* IGEQ * KGPRNT) 

GO TO 110 

CALL PATHS(ICALC»NERRQR*IGEOjp KGPRNT) 

IF (NERR0R.GE.2 ) GO TO 160 
IF U.LT.WPL) GO„TD 130 
CONTINUE 

-COMPUTE .JJSED* RESERVE* TQTAL* .AND INTERNAL FUEL WEIGHT. 
UFUSED-WGTD-H 

IF (IhS.EQ.O) WFUSED-WFUSED-WM1SS 
IF (IBS.EO.O) WFUSED^WFUSED— WBOMBS 
IF ( lAS.fcQ.O) WFUSED*WFUSED— WAMMUN 
IF (ITS.EQ.C) WFUSED-WFUSED-WETANK 
WFUSED a WKFUEL*WFUSED 
WFRES-WFUSED/ (l.-FRFURE)-WFUSED 
WFTLT«WFUS£C+WFRES+WFTRAP 
WFUEL-WFTOT-WFEXT 
IF (W.GE.WPLJ GO TO 140 
IF (KERRUR.E0.2) WRITE (6*602) 

FORMAT ( /37H FUEL WEIGHT EXCEEDS AVAILABLE HEIGHT/) 
WFUtL-WGTO-HPL 
GO TCJ 160 

-CALCULATE LANDING FIELD LENGTH. 

RHO". 00238 
I A(i M 12 

CALL TRAJ0K2) 

CALL AERO 1 1 CALC* NERROR* IGEQ* KGP RNT ) 

IF (NERRQR.GE #2) GO TO 160 
CALL TRAJCK1) 

WLAND«WGTQ-WKLAND*WFUSED 
VSTALL“SQRT(2.*HLAND/(SWING*CL*RH0) ) 

VSCREN-1.3+VSTALL 
VTD0WN*1.15*VSTALL 

XAIR« ( [ VSCREN**2-VTD0WN*+2) / 64,4+50. )*RLD 
XGRLAN*VTOOWN++2/(64.4*DECEL) 

FLLAND»(XAIRfXGRLAN)/FLFAC 
-CALCULATE ENDURANCE ALTITUDE AND TIME. 

SHN-HENDUR 

CALL GOLDEN (2000.* 55000. *2* 12* XX* Y* NERRQR* TGEO* KGPRNT J 
IF (NERRQR.GE. 2) GO TO 160 
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116 HENDUR-XX 

TtNDUR*Y+ALOG(WGTO/{WGT O-WFUSED ) ) 

C CUMPUTE GROSS WEIGHT. 

WGCALC*WFIXED + WPl.+ WFTOT 
IF ( IPRINT.EQ.O) GO TO 16G 

1ZC WRITE ( 6# 603 \ WGTO# WFUSED# WFRES# WFTOT, WFUEL# WPL# W» 

1 wLAND#XGRLAN#FLLAND# HENDUR# TENDUR# WGCALC 
603 FCRMATl///7X#7HLANDING///3X#6HWGT0 «E14. 7# IX# 8HWFUSED -E14.7# 

1 2X#7HWFRE5 “E14»7#2X#7rt WFTOT «E14.7#ZX#7HWFUEL -E14.7/ 

2 4X# 5HWPL «£14.7»6X#3HW -E14. 7# ZX, 7Hw(.AND «tl4.7# 

125 3 IX# 8HXGRL AN * E14. 7# IX# 8HFLL AND *E14 .7 /IX# 8HHENDUR -E14.7# 

4 IX# 8HTENDUR *E 14. 7# IX# 8HWGCALC -E14.7/1 
160 RETURN 

END 
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FTN 4.5+410 05 


v 



1 


5 


10 


15 


2C 


25 


30 


35 


4C 


45 


5C 


55 


SUBROUTINE TAKEOF ( ICALC, NERROR, IGEO, KGPRNT ) 

COMMON /TRAJCH/ ALPHA, ARW, CD, CD L, CDO,C L,DESLF, DR AG, EN,HN, KP, 

1 ft LD, SMN, QMAX, RANG E,SFC, SWING, THRU ST »Tw,ULTLF,W,WAMMUN, WE TANK, 

Z wFUEL, WGTOWT,WMISS,WTOT,WPLWT,HSTART(1Z),HSTART(12) ,TIH( 12) , 

3 PSlGT(xZ),TDOTSH12>,NZST(12),PSIT(12),TDOTIT(12),NZIT(12), 

4 W BOMBS, wFEXI,SFCU, THRSTU, GDI NS P, PR TOT, WKFUEL,CRHACH, FLTO, XI 12), 

5 FLLAND,TENDUR, BLR ANG, BLT IME, TI MT01, T I MT02, WFTO, PNAME1 ( 12 ) , 

6 PNAME2(12),SHNT(12),HNT(1Z),WFT(12)>TIHET(1Z),XT (12),CLT(12), 

7 CDT(12),ALPHAT(12>»CLIT(12),CDIT(1Z)» ALPHIT 112 } ,WFUSED, WLAND, 

8 QT(12) >ftLDT(12>,$FCT(12),TNT(lZ),CET( 12),DY(18), 

9 IAO,IPS,ITS,IWS,IBS,NPHASE, IDY(14) 

COMMON /TRAJEX/ ALP1GT l 12 ) , C DINST ( 12) , CDIGTt 12 ) , CL1GT ( 12) , DECEL, 

1 DELHP»DfcLMP,DXCRUS, FLFAC, FRFURE, HEND ( 12) , HENDUR ,HHAXP ,HMINP, 

2 HPATHT ( 50), MEND (12) , MENDUR,NT( 12) , PRTOTTt 12 WRADIT(IZ), 

3 RADST (12) ,SF,SFCUT(12), SHMAXP, SMHINP, 

4 TNUT(12),VIN0(1Z),VELT(12),VPATHT( 50) , WARYC 12 ) , WCOP.BP, WEELT (12) , 

5 WFRES,WFTOT,WFT01,WFT02,WFTRAP,WKLAND,XUESC,XGRLAN, FWGMAX,TOL, 

6 WFI XEO , WGC AL C> WP L, WGT □, DUMY { 8 2 ) , I A ( 12 I, I/iS, IB (12 ), I BREG, IENDUR, 

7 IMISS, IP ( 12) , IPHASE, IPLOT, 1PRINT, I PRT ( 12 ) » I PSIZE, IP STQ1, IPST02, 

8 IW(12),IX(IZ),KERRGR,MILC0M,NCRU5E,NLEGCL,NLEGCR»NLEGL0, NHISS, 

9 NPATH, IDUMY ( 6 ) 

DIMENSION ANS ( 4 ) 

C —CALL PROP TO SIZE ENGINE. 

IPS*IP$IZE 
KN-O. 

SMN-U. 

CALL TRAJ0K2) 

CALL PROP ( I CALC, NERROR, IGEO, KGPRNT) 

IF (NERR0R.GE.2) RETURN 
CALL TRAJOl(l) 

C CALCULATE STARTUP AND TAXI FUEL USED. 

IPS-IPST01 
CALL TRAJ0K2) 

CALL PROP (1CALC, NERROR, IGEO, KGPRNT) 

IF (NERR0R.GE.2) RETURN 
CALL TRAJOl(l) 

TN-THRUST+EN 

WFTQI*. 0l67*TIHTbl*SF'C*TN 

C BEGIN FIELD LENGTH CALCULATION. 

HNTO-HSTART(l) 

HN-HNTO 

CALL AT62(HN, ANS) 

RHO-ANS(l) 

A*ANS<4) ' 

C GUESS STALLING HACH NUMBER. 

SMN-.25 
IAO-14 

CALL TRAJOl(H) 

CALL AEROdCALC, NERROR, IGEO, KGPRNT) 

If (NERROR. GE. 2) RETURN 
CALL TRAJOl(l) j 

CLS-CL i 

C CALCULATE STALLING SPEED. T 

\1S«SQRT(2.*WGT0/(RH0*CLS*SWING) ) j 

V2*=l.ms 

SMN2-V2/A 
SMN-SMN2 
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115 


CALL TRAJ0112) 

CALL AERO(ICALC#NERROR# IG EO# KGPRNT) 

IF ( NERROR. GE. 2) RETURN 
CALL TRAJ0H1) 

CL2*CL 

C FIND THRUST AT END OF TA (EQFF. 

IPS-IPST02 
CALL TRAJD112) 

CALL PRDP( 1C ALC# NERROR# IGEO#KGPRNT) 

IF ( NERROR . GE . 2 ) RETURN 
CALL TRAJOl(l) 

TN*THRUST*EN 

TN2«TN 

SFCZ-SFC 

C FIND STATIC THRUST. 

SHN-O. 

CALL TRAJ0K2) 

CALL_ PROP (ICALC# NERROR# IGEO# KGPRNT) 

IF (NERROR. GE. 2) RETURN^ 

CALL TRAJOl(l) 

TN«THRUST*EN 

TNO-TN 

SFCO-SFC 

C -CALCULATE AVERAGE THRUST AND SFC DURING TAKEOFF. 

TNAVE«.5*t TN0+TN2) _ 

SFCAVE-.5* (SFC0+SFC2) 

C CALCULATE TAKEOFF FIELD LENGTH REQUIRED. 

IF (MILCOM.EQ.O ) TOOBHT-5D.O 
IF (MILCOM.NE.Ol TOOBHT-35.0 

FLTO«2.1U*(.37^*T006HT+.01163*MGTO/(RH0*CL2*SHINGn* 

1 (2.7+l.O/lTNAVE/WGTO— .04)) + 32.0/,SQRT( RHO ) 

C CALCULATE TAKEOFF FUEL USED. 

IF (TIMTO2.LT.0.Q) GO TO 10 
IF CTIHT02.EQ.0.0) GO TO 20 
HN-O. 

IPS-IPSTQ2 
CALL TRAJ01C2) 

CALL PROP (I CALC#NERROR#IGEO# KGPRNT) 

IF (NERROR. Gt. 2) RETURN 
CALL TRAJOl(l) 

TN-THRUST+LN 

WFT02».0167*TIMTQ2*$FC*TN 

GO TO 30 ... . 

10 WFT02»TNAVL+SFC*FLT0/ (1800. *V2) 

GO TO 30 
20 WFT02*0.0 

C CALCULATE TOTAL FUEL USED FOR STARTUP# TAXI# AND TAKEOFF. 

30 VFT0-HFT01+RFT02 

W*WGT0-WFT0 

IF (IPRINT.EQ.Q) RETURN 

WRITE ( 6# 60C ) WGT0#WFTQ1»WFT02#WFT0#W#HNT0»CLS# VS# V2# SHN2#CL2# 

1 TN2#SFC2» TNO#SFCO#TNAVE#SFCAVE#FLTO 
600 FQRHAT(1H1#6X#7HTAKEOFF///3X>6HWGTO «£14.7#2X# 7HWFT01 -E14.7# 

1 2X#7HWFT02 «E L A . 7# 3X# 6HWFTD «E 14. 7# 6X» 3HW -E14.7/ 

2 3X# 6HHNT0 - E14 . 7» 4X# 5HC LS -E14. 7# 5X#4HVS -E14.7# 5X# 4HV2 -E14.7# 

3 3X#6HSMN2 *E14 . 7/ 4X# 5HCL2 -E14. 7# 4X#5HTN2 -E14.7# 

4 3X» 6HSFC2 m fc 14.7# 4X# 5HTN0 -E14. 7# 3X# 6HSFC0 *E14.7/ 

5 2X# 7HTNAVE -E14.7# IX# 0HSFCAVE -E14. 7# 3X#6HFLT0 -E14.7/) 

RETURN 

END 
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SUBROUTINE CLIMB(1QALC*NERR0R*IGE0*KGPRNT) 

REAL MSTART*MEND 

COHMON /TRAJCM/ ALPHA* ARW*CD*CDL»CDO* CLjsDESLF*DR AG* EN*HN* KP* 

1 RLO* SMN*QMAX> RANGE* SFC* SWING* THRUST*TW*ULTLF*W*WAHHUN*WETANK* 

2 WFUEL* WGT0WT*WMISS*WTDT*WPLWT»MSTART(12I* HSTART (12)*TIM(12)* 

3 PS1GT(12)*TD0TST(12)*NZST(12)* PS IT (12)* TDQT IT (12)* NZIT( 12 )» 

4 WBOMBS»WFEXT*S p CU*THRSTU*CDINSP*PR TOT * WKFUE L * CRMACH *F LTQ* X( 12) > 

5 F LL AND * TENDUR » BL RANG* BLTIHE*TIHTD1*TI MTQ2# WFTO* PNAHE1 (12 ) » 

6 PNAME2(12)*SMNT(12)*HNTUZ)*WFT(12)*TIMET(12>*XT(12)*CLTU2)» 

7 CDT(12)*ALPHATilZ)*CLIT(12)*CDIT(12)pALPHIT(12)»WFUSED*WLAND* 

B QT(12)*RLD7(12)*SFCT(12)*TNT(12)*CET(12)*DY(18)» 

9 IA0»IPS*ITS»IWS*IBS*NPHASE»IDY(14) 

COMMDN /TRAJEX/ ALPlGTU2K*gDINSTU2)*C01Gm2)»CLlGT(12)*DECEL* 

1 DELHP*DELMP*DXCRUS*FLFAC*FRFURE*HEND(12I*HENDUR* HMAXP*HHINP* 

2 HPATHT(501*HENDU2)*HENDUR*NT(12)*PRT0TT(1Z)*RADJT(12)* 

3 RADST(12)*SF*SFCUT(12)*SMHAXP*SMHINP* 

4 TNUT(12)*VIND(12)*VELT(12)*VPATHT(50)*WART(12)*WC0HBP*WEELT(12)* 

5 WFRES* WFTDT* WFT01* WFT02* WFTRAP* WKIAND * XDESC * XGRLAN*FWGHAX*TQL* 

6 WFIXED*UGCALC*WPL*WGT0*DUNY(BZ)*IA(1Z)*IAS*IB(12)*IBREG*IENDUR* 

7 IMIS5*tPJ12)*IPHASE*IPLDT*IPRINT»IPRT(12)*IPSIZE»IPST01*IPSTQ2* 

B IW(1Z)*IX(12)*KERRQR*MILC0M* NC RUSE* ML EGCL* N LEGO R*NLEGLO* NKISS* 

9 NPATH* IDUMYC6) 

DIMENSION ANS ( 4 ) 

IF ( IPRT(IPHASE).NE.Q) WRITE (6*600) 

600 FORMAT ( / / / 18H CLIMB PHASE//) 

C FIX STARTING ALTITUDE 

IF (HSTARIilPHASE) .GE.O.O) HNINT"HSTART( I PHASE) 

IF (HSTART(IPHASE).LT.O.O) HNINT-HNT (I PHASE- 1 ) 

C FIX ENDING ALTITUDE 

HNFIN-HEND(IPHASE) 

IF ((IBREG.EQ.0).OR.(HEND(IPHASE).GT,0.0)) GO TO 40 
SHN-MENDCIPHASE) 

CALL GQLDEN(2000.* 60000, *1*12* XX* Y* NERROR* IGEO*KGPRNT) 

IF (NERROR.GE.Z) RETURN 
HNFIN*XX 

4C SMN-MSTARTUPHASE) 

IF (HSTART(IPHASE) .LT.O.O) SMN-SMNT (IPHASE-1 ) 

IF (HNFIN.LT. 80000.1 GO TO 43 

NLEGS-20 

GO TO 45 

43 NLEGS«20.*ALQG( ( 4000.-.05+HNINT ) / (4000 .-.05+HNFIN ) ) + . 5 
IF (NLEGS.LT.3) NLEGS*3 
IF (NLEGS.GT.20) NLEGS-20 
45 IF (NLEGCL.GT.O) NLEGS-NLEGCL 

NLEGS1-NLEGS+1 
B1«HNFIN-HNINT 
B2«1.718262/FLOAT(NLEGS) 

IF (Bl.GT.L.O) GO TO 50 
SMNTt IPHASE ) *SMN 
HNTCIPHASE J-HNFIN 
CLT( IPHASE )*Q.O 
ALPHAT(IPHASE)«0.0 
NFT ( IPHASE ) *C *0 
TIMET(IPHASE)»0.0 
SFCT ( IPHASE ) ■ 0» 0 
TNT ( IPHASE ) «0 » 0 
CDT ( IPHASE ) *0. Q 
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RLDT{IPHASE)-0.0 
WtELT( IPHASE ) -W 
WART(IPHASE)*0.0 
RETURN 

50 HN-HNINT 

TIME-0.0 
WF*0.0 
XCLIMB-0.0 
VINDFS-0.0 
IPS-IP ( IPHASE ) 

IAO-8 

IF (VIND(IPHASE).IE.O.) GO TO 60 

C INITIALIZE CONDITIONS FDR CONSTANT INDICATED AIRSPEED CLIMB. 

IF ( IPRT(IPHASE) .NE.O) WRITE (6,601) 

601 FCRMATI/34H CONSTANT INDICATED AIRSPEED CLIMB//) 

IF (MSTART(IPHASE) .GE.O.) VINDFS*1.609*VIND( IPHASE J 
IF CMSTART(IPHASE) .LT.O. ) VINDF5«VSLT( IPHASE-1) 

CALL AT 62( HN, ANS ) 

RHOP-ANS (1 ) 

CL«2.*W/(.002378*VINDFS**2*SWING) 

CALL TRAJ0112) 

CALL AEROC 1C ALC, NERRQR, IGE0,KGP RN T ) 

IF ( N ERROR * GE .2 ) RETURN 
CALL TRAJOl(l) 

CALL TRAJ0K2) 

CALI PROP{ICALC,NERROR,IGEO,KGPRNT> 

IF (NERR0R.GE.2) RETURN 
CALL TRAJOl(l) 

TN-THRUST *EN 
TNU*THRSTU*EN 

60 IF (IPRT(IPHASE).EQ.O) GO TO 70 

WRITE (6,602) W,HNINT,HNFIN, SMN, HN 

602 FORMAT £ / 6X,3HW *E14.7,2X,7HHNINT »E14.7,ZX,7HHNFIN -E14.7, 

1 4X,5HSMN -E 14. 7»5X, 4HHN -E14.7/) 

70 DO 140 LEG-1, NLEGS1 

IF (LEG.LT.NLEGS1) "GO TO 80 

IF {IBREG,EJ3.0.QR.HEND{IPHASE).GT.0. ) GO TO 140 

HN1-HN 

CL1-CL 

ALPHA1-ALPHA 

SFC1-SFC 

TNI* TN 

CD1-CD 

WTOT1-WTOT 

HL0W-HN1 

HHIGH-HN1+3000. 

SMNSAV-SMN 

SMN- AM INI (SMN, MEND (I PHASE) ) 

CALL GOLDEN (HLOW, HHIGH, 1, 6, XX, Y, NERRDR, IGEO, KGPRNT) 

SMN-SMNSAV 

HNFIN-XX 

DELH-HNF IN-HN1 

HN-HN1 

CL-CL1 

ALPHA-ALPHAl 
SFC-SFC1 " 


ORIGINAL PAGE IS 
OF POOR QUALITYj 


SUBROUTINE CLIMB 


76/76 OPT-2 
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CD-CD1 

WTOT-WTOT1 

IF CNERRQR.LT. 2. AND. DELH. GT. 0. 005+HN ) GO TD 81 

NERROR-0 

GD TO 140 

80 CiELH-Bl+ALOG( C 1 . +FLQAT ( LEG J *B2 3 / C 1 . +FL OAT t LE G— 1 1 *B2 ) ) 

81 IF < VINDFS.GT .0 . ) GO TO 120 
CALL AT 62 ( HN> ANS > 

RH0-ANSC1) 

A* ANS (4) 

IF (LEG.EQ.l) VEL-SMN+A 
IF ( LEG . GT ■ 1 ) S MN* VEL /A 
IF CSMN.LT .0.0) NERROR-2 
IF CSMN.LT. C.O) WRITE (6*610) IPHASE* LEG 
610 FORMAT C/////34H *****FATAL ERROR IN CLIMB, PHASE«*I2*7H 

1 27H **+**NEGATIVE MACH NUMBER./////) 

IF CSHN.LT.C.O) RETURN 
Q-.5*RH0*tfEL**Z 

QS-Q+SWING ... 

CALL TRAJ31C2) 

CALL PROPl ICALC*NERROR* IGEO*KGPRNT) 

IF (NERROR.GE.Z) RETURN 
CALL TRAJ01C1) 

TN-THRUST+EN 

TNU«THR$TU*EN 

IF (TN/W.LT.1.5 ) GO TO 90 

C EXECUTE FOLLOWING SECTION ONLY IF TN/W ,GE. 1,5 

SING-1.0 

IA0-0 

CALL TRAJ01C2) 

CALL AEROCICALC>NERROR>IGEO>KGPRNT) 

IF CNERR0R.GE.2) RETURN 
CALL TRAJOltl) 

CD-CDO 
GO TO 100 

C EXECUTE FOLLOWING SECTION ONLY IF TN/W .LT. 1.5 

90 CL*U/aS)+SQRT( 1.-. 444* C TN/W) **2> 

IAO-8 

CALL TRAJ01C2) 

CALL AEROCICALC»NERRQR*IGEO»KGPRNTJ 
IF CNERROR.Gt .2 ) RETURN 
CALL TRAJOltl 1 
CL-CL-TN+SINC ALPHA/57.3) /QS 
CALL TRAJJ1C2) 

CALL AEROCICALC*NERROR>IGEO,KGPRNT) 

IF (NERROR.GE.Z) RETURN 
CALL TRAJOl(l) 

EMAX-. 5+SQRT (CL**2/(CD0*CCD— CDU) ) ) 

Al-TN+EMAX/W 
A2-S0RT ( Al**2 + 3 * ) 

SING-.B333*TN/W— .167*A2/EMAX— 1.5/CEMAX+CA1+A2)) 
IFCSING.LE.O.O) SING-0.0 
1FCSING.LE.O.O) GO TO 110 
COSG-SQRTC 1 .0— SING* *2 ) 

ICO DELT-DELH/(VEL*SING) 

IF (DLLT.GT.0.0 ) GO TO 110 
NERROR-2 
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0 


175 


L80 


L85 


L90 


195 


*00 


'.06 


*10 


:is 


•2C 


:Z5 


teR IT£ (6* 602 ) IPHASE# LEG 

6Q3 FQRHAT i /// / /34H ***#*FATAL ERROR IN CLIMB. PHASE«#I2#7H LEG-#IZ/ 
1 26H *****NEGATIVE DELTA TIME./////) 

RETURN 

110 DVDT«32.2*nN-CD*QS-W*SING)/W 

IF tDVDT.LT.O.) SING- AMAXlt 0. # ( TN-CD*QS ) /W ) 

IF tSINS.LE.O.) DVDT-32.2*(TN-CD*QS)/W 
DRaG#CD*QS 

IFtSING.LE.O.O) DELT-60.0 
DELV«DVDT*DELT 
IFtSING.LE.O.O) C0SG-1.0 
IFtSING.LE.O.O) DELH*0 »Q 
DELX»VEL*DELT*CQSG/6Q80. 

VEL-VEL+QELV 

DELWF«TN*SFC$DELT/3600. 

C TIME TO CLIMB 

TIME«TIME+DELT/60. 

C FUEL USED DURING CLIMB 

WF*WF+DELWF 

C GROUND DISTANCE COVERED DURING CLIMB 

XCLIMB-XCLIMB+DELX 

HN-HN+DELH 

W*W-DELWF 

IF (W.GT.WPL) GO TO 115 


IF tKERRQR.IQ.2J MRJtTE_ t6±605) IP BASE# LEG_ 

605 FORMAT t / /26H W.lcoWPL IN CLIMB. PHASE-#I2#7H LEG## 12/) 

RETURN . _ . 

115 IF tlPRT(IPHASEi.EQ.O) GO TO 1.40 

WRITE {6#606)_LEG#SHN#HN#RH0#a#VEL#SFC#TN#ALPHA#CL»CDD#CD#EHAX#A1# 
1 A2, SING# CQSG#DELT#DVDT#DELV# DELX#OELWF# W»WF# TIME# XCLIMB 

606 F DRHATt /4X#5H LEG «I2/4 X#5HSMN «E14.?#5X#4JHN »E14.7#4X» 

1 5HRHD »E14.'f#6X#3HA *E14.7,4X# 5HVEL «E14.7/4X#5HSFC »E14.7#5X# 


2 4HTN #E14.7»2X»7HALPMA_#E14*7»5X#4HCL «E14.7»4X#5HCDO #E14.7/5X# 

3 4HCD «E14.7#3X#6HEMAX #E14. 7#5X# 4HA1 #Ei4.7# 5X#4HA2 #Ei4.7#3X, 

4 6HSING «E14.7/3X#6KCDSG ^14.7# 3X# 6HDELT #E14.7#3X#6HDVDT -E14.7# 

5 3X# 6HDELV «t 14 .7# 3X# 6HDELX »E14.7/2X#7HDELWF »E14.7#6X#3HW #» 

6 E 14.7# 5X #4HuF «E14.7»3X#6HT IHE »E14.7#2X#7H XCL IM B»£1 4.7/ ) 

GO TO 140 

C " START OF CONSTANT INDICATED AIRSPEED SECTION 

120 hn-hn+deLh ~ 

CALL AT 62( HN# ANS ) 

RHQ#ANS 1 1 ) 

A#ANS 1 4 ) 

Cl#. 5* t RHO+RHDP) 

C2*SQRT tCl ) _ 

C3*20.51*C2~ 

C4-C1+C2 

C5-32. 2*C3+D£LH/VINDFS-.02438*VINDFS*l RHQ-R hQP) /C4 
C6«32.2*(TN/W-.0011B9*\/INDFS*+2*CD*SWING/W) 

DELT-C5/C6 

IF lDELT.GTiD.OJ GQ_TQ L25 _ __ 

NERROR# 2 

WRITE (6# 607) IPHASE#LEG#DELH#V1NDFS#RH0#RH0P#TN#W#CD#SWING#DELT# 

1 CI#C2#C3»C4#C5#C6 

607 FDRMAT {/Z///34H **+**FATAL ERROR IN CLIHB. PHASE-#I2#7H LEG-,12/ 
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0 


230 


235 


240 


245 


25C 


155 


260 


265 


270 


275 


28C 


285 


125 


130 


606 


C — 

c 

c — 

14C 


c 


1 40H *****NEGATIVE DELTA TIHE IN VIND CLIM3. /5X* 5HDELH=,E12. 5* 2X* 
Z ?HVlNDFS**E12.5»2X*4HRHa«,E12.5/5X*5HRHaP«*E12.5»2X*3HTN«,E12.5* 

3 2X>2HW«,E12.5/5X*3HCD*»E12,5>2X,6HSWING»*E12.5*2X,5HDELT*,E12.5/ 

4 5X* 5HC 1-C6/E13* 5/ E13.5/E13.5/E13.5/EI3.5/E13.5/////) 

RETURN 

TIME-TIHE+DELT/60. 

VEL*VINDFS/SQRT(RHO/. 002378) 

SMN=VEL/A 


IF (SHN.LT.O.O) NERR0R«2 
IF tSMN.LT.O.O) WRITE (6*610) I PHAS E*L EG 
IF (SMN.LT .0.0) RETURN 
SING«DELH/ ( VEL*DELT) 

IF (SING.GT.l. ) 5ING«1.0 ... 

C0SG«SQRT(1.-SING**2) 
DELX*VEL't ! 0tLT*COSG/6O8O._ _ 
XCLIHB-XCLIHB+OELX 

CALL TRAJ01U) 

CALL PROP(ICALC*NERRQR*IGEO*KGPP,NT) 

IF t NERROR. GE « 2) RETURN 

CALL TRAJ01 (1 ) 

TN-THRUST*EN 


ORIGINAL PAGE IS 
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TNU*THR5TU*EN 
DELWF»TN*SFC+DELT/ 3600. 


WF*WF+DELWF 

W«W“DELWF „ 

IF (W.GT.WPL) GO TO 130 

IF (KERR0R.EQ.2) WRITE (6*605) IPHASE*L£G 
RETURN 

C 7*1 •-( DEL H/ ( VINO FS+DELT+C3 ) ) **2 
IF (C7.LT.G.) C7-0.0 

CL«H.*W*SQRT (C7) /( .002378*VINDFS**2*SWING) 

IAQ*8 

CALL TRAJ0K2) 

CALL AERO! ICALC*NEKROR, IGEO*KGPRNT) 

IF (NERROR »GE. 2) RETURN 
CALL TRAJOl(l) 

RHQP-RHO 

IF (IPRT(IPHASE)eEO.O) GO TO 140 

WRITE (6,6GB) LEG* VINDFS* SMNi-HN* DELT, DELX*C6* VEL,DELWF *W»TN*CL*CD 
FORMAT(/4X,5HLEG "12/ 

1 IX* 8HV INDFS «E14.7*4X,5HSHN «E14.7* 5X,4HHN -E14.7, 

2 3X*6HDELT -E 14. 7, 3X*6HDELX *El*t. 7/ 5X* 4HC6 -E14.7* 

3 4X> 5HVEL “E14. 7* 2X, 7HDELWF «S14. 7, 6X, 3HW -E14.7* 

4 5X*4HTN ■ E14« 7/ 5X* 4HCL » E14. 7, 5X* 4HC0 «E14.7/> 


END 0? CONSTANT INDICATED AIRSPEED SECTION 


CONTINUE 

CALL AT62( HN* ANS ) 

A« ANS ( 4 ) 

SMN-VEL/A 
RHQ*ANS (1) 

Q"«5+RHQ*VEL*42 

DXCRUS*OXCRUS + XCL IHB *FLOAT(IXtIPHASfc ) ) 
BLTIME"BLrif»E + TIHE/60. 
BLRANG*8L«ANG+XCLIMB 
TEST FOR WEAPONS DROP 
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05 


IF (IWdPHASE) .EO.O) GO TO 150 

W-W-fcHISS 

IWS-0 

C TEST FOR BOMBS DROP 

!90 15C If ( IB 1 IPHASE ) « EO.O ) GO TO 160 

Vt«k-WB0HBS 

IBS-0 

C TEST FOT AMMG DROP 

160 IF ( I A 1 1PHASE ) « EQ.C I GO TO 170 
295 W“W— WAMMUN 

I AS-G 

170 IF ( IPRTl IPhASE) .EO.O ) GO TO 1B0 

WRITE [ 6j 609 1 SHN*HN*VELj’XCLIMti>DXCRUS»TlME#WP*W 
609 FORMAT t /4X, 5HSKN *E14 .7, 5X, 4HHN =£14. It 4X* 5HVEL “E1A.7, 
30C 1 1X.BHXCUNB =£ 1 A. 7, 1 X» BHD XCRUS -El A. 7/3X ,6H T IME =E1A,7> 

2 5X*4HWF *E1A.7»6X>3HW -E14.7/) 

160 SHNT(IPHAS£)«5MN 
HNTt IPHASE) =HN 
CLT( IPHASE) «CL 

305 ALPHATt IPHASE '-ALPHA 

WFTUPHASE )-Wf 
TIMETUPHASE J-TIHE 
VELT l IPHASE ) -VEL 
SFCTIIPHASE)»SFC 

310 TNT t IPHASE ) -TN 

CDT l IPHASE J-CO 
RLOTIIPHASEJ-RLD 
REEL T( IPHASE 1 »H 
WARTtlPHASE)*WTQT 
315 QT ( I PHAS E ) » Q 

SFCUTIIPHASEl-SFCU 
7NUTIIPHASEJ-TNU 
CDINSUIPHASEJ-CDINSP 
PRTOTTt IPHASEJ-PRTOT 
3Lu XTUPHASEJ-XCLIMB 

RETURN 

END 


■: 1 

;! 


y 

I 




SUBROUTINE CRUISE 


76/76 OPT-2 


FTN ‘t. 5+410 


1 


5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


55 


600 
C— - 


1C 


SUBROUTINE CRUISE{ICALC,NERRDR» IGEO,KGPRNT) 

REAL HSTART 

COMHON /TRAJCM/ ALPHA, ARW, CD, CDL,CDQ, C L, D ESL F, DR AG, EN, HN, K P, 

1 RLD,SNN,OftAX, RANGE, SFC, SUING, THRUST,TW,UL TLF, W, WAHMUN, WE TANK, - 

2 WFUEL , WGTOWT , WHISS, WTOT,WPL WT ,HSTART(12),HSTART(12),TIM(12) » 

3 PSlGTll2},TD0TST(12),NZST(12),PSir{12),TD0TIT(i2J,NZIT(12), 

4 WBOMBS,WFEXT,SFCU,THRSTU,CDINSP,PRTOT»WKFUEL,CRMACH,FLTD,X(12), 

5 FLL AND, TENDUR, BLRANG, BLTIME,TIMT01,TIMTD2,WFTO,PNAME1(12), 

6 PNAME2 ( 12) ,SMNT(12) ,HNT (12J,WFT(12),TIMET(12),XT(12),CLr(121, 

7 CDT{12),ALPHAT(12),CLITC12),CDim2),ALPHIT(12),WFUSED,WLAND, 

8 QT(12),RLDT(12),SFCT(12), TNT (12),CET(12),DY(18), 

9 I AO, IPS,ITS,IHS,IBS,NPHASE,IDY<14> 

CDHHON /TRAJEX/ ALP1GTI 12 )# CDINST ( 1Z ) , CD1GT ( 12 ) , CL 1GT < 12) , OE CEL, 

1 DELHP,DELMP,DXCRUS,FLFAC,FRFURE,HEND{ 12) , HENDUR , HMAXP,HMINP, 

2 HPATHn50>,HENDC12>,MENDUR,NT{12),PRT07TC12),RADIT(lZ), 

3 RADST(12),SF,SFCUT(12)»SHMAXP, SMHINP, 

4 TNUT(12),VIND(12),VELT{12),VPATHT(50),WART(12),WC0HBP,WEELT(12)f 

5 WFRES,WFTOT,UFT01,WFTaZ,WFTRAP, HKL AND, XDESC, XGR LAN, FWGMA X,T£JL, 

6 WFIXED,wGCALC,UPL,WGT0,DUHY(82 ) , I A{ 12 ) , I AS, IB( 12) »I BR EG, 1ENDUR, 

7 IMISS,IP(12) , I PHASE, I PLOT »IPRINT,IPRT(12),IPSIZE,IPSTD1,IP$TD2, 

8 IW(I2),IX(12),KERR0R,HILC0K,NCRUSE,NLEGCL,NLEGCR,NLEGL0,NHISS, 

9 NPATH, IDUMY(6) 


DIMENSION A NS ( 4) ,HNTAB{ 10) ,SHTAB(10) 

DATA HNTAB/0.,5000.,10a00., 15000., 20000., 25000., 30000., 

1 35000. ,37000. ,100000./ 

DATA SHTAB/34500. 732900, ,31400. ,30500., 29400., 28 ZOO, ,27110.* 
1 25950, ,20900. ,20900./ 

DATA D MON/O.O/ 

IF (IPRT(IPHASE).NE.O) WRITE (6,600) 

FORMAT ( ///19H CRUISE PHASE//) 

—TEST FDR TANKS DROP_ 

WFUSED-WGTD-W 

IF (IWS.EQ.O) UFUSED-HFUSED-WMISS 
IF (IBS.EQ.O) WFUSED-WFUSED-WB0M8S 
IF (IAS.EQ.O) WFUSED-WFUSED-UAMMUN 
IF (ITS.EQ.l.AND.WFUSED.GE.WFEXT) W-W-WETANK 
IF ( ITS • EQ. 1. AND. UFUSED. GE . WFEXT ) ITS«0 
IPS*IP ( IPHASE ) 

SMN-MSTART(IPHASE) 

HN«HSTART( IPHASE) 

IF (MSTART(IPHASE).LT.O.O) SMN*SHNT ( I PHASE-1 ) 

IF (HSTARr(IPHASE).LT.O.O) HN«HNT ( I PHASE-1 ) 

IF t IBREG.EQ. O.DR.HSTAP<T (IPHASE).NE.Q.O) GD TQ 


obigwai Mj® ® 

OF POOR QUAIJ-W- 


10 


CALL GDLDEN(2006.i60dd0.,l,12,XX,Y,NERR0R;iGE0,KGPRNT) 
IF (NERRDR.GE.2) RETURN 

HN«XX 

CALL AT62(HN,ANS) 

RHD»ANS(1) 

A»ANS ( 4 ) 

VEL* SMN* A 

CL«2.*W/(RHO*SWING*VEL**2) 

IA0*8 

CALL TRAJ0K2) 

CALL AERD( IC ALC , NERROR, IGEQ,KGPRNT) 

IF (NERRDR.GE.2) RETURN 
CALL TRAJOl(l) 

CQSA*CUS(ALPHA/57.3) 
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05 


6C 


65 


7C 


75 


60 


05 


9C 


95 


IOC 


105 


110 


TIME-G.O 

WF-0.0 

IF ( X ( IPHASE) .GE.O.)XCRUSE»X( IPHASE)+OXCRUS+XDESC*FLOAT (IX ( IPHASE) ) 
IP ( Xt IPHASE) . LT ,0. ) XC RUSE “RANGE /FLOAT ( NCRUSE 1 +DXCRUS+ 

1 XOESC+FLDAT (IX(IPHASE) ) 

IF (XCRUSE ,LE. 0.) DXCRUS-Q.O 
IF (XCRUSE .LE. 0.) GO TO 90 
IF ( X C I PHASE; 1 .LT.100. > NLEGS-2 

IF ( (X(IPHASE) . GE • 100 « ) .AND* £ X ( IPHASE) .IT • 30 0 • J ) NLEGS-3 
IF ((X(IPHASE). GE. 300. KANO. (X(IPHASE) ,LT. 1000, )) NLEGS-4 
IF ((X(IPHASE).GE. 1000.). AND. (X(IPHASE).LT. 3000.)) NLEGS* 5 
IF (X(IPHASE) .GE.3000. ) NLEGS-6 
IF (NLEGCR.GT.O) NLEGS-NLEGCR 
DELX-XCRUSE/FLOAT( NLEGS) 

IF (IPRT(IPHASE).EQ.O) GO TO 20 

WRITE (6# 601) W.SMN.HN.RHO, A, VEL.CL, ALPHA# TIME.WF, XCRUSE. DXCRUS, 

1 DELX » I AO, IPS 

601 F CRMAT ( /6X, 3HW -E14 . 7, 4X, 5HSMN -E14.7, 5X, 4HHN -E14.7,4X« ■ «H0 * 

1 E14. 7. 6X» 3HA »E14 . 7/4X, 5HVEL - E14. 7, 5 X, 4HCL -E14.7.2X.7 ._PHA - 

2 E14.7, 3X.6HTIME -E14. 7, 5X, 4HWF -E14 . 7/ IX. 6HXCRUSE -E14.7.1X, 

3 BHDXCRUS -E14.7.3X.6HDELX -E14.7.4X.5HIAQ -12, 16X, 5HIPS -12/) 

20 DO 50 LEG-1, NLEGS 

IF (HEND(IPHASE).EQ.-l.) GO TO 21 
RHOl-RHO 

RH02-2.*W/(CL*SWING*VEL#*2) 

CALL TAINT (HNTAB, SHTAB, HN, SB# 10, i.NERR.OMON) 

HN-HN+2.+SH* (RH01— RH02 )/ (RH01+RHQ2) 

CALL AT62(HN, ANS ) 

A-ANSI4) 

VEL-SHN*A 

RH0-2.*W/(CL+SWING*VEL**2) 

C«.5*RH0*VEL**2 
QS-Q+SWING 
GO TO 22 

21 HN«HSTART( IPHASE) 

CALL AT62(HN,ANS) 

A-ANS14) 

KHO-ANS (1) 

VEL-SMN*A 
G-.5*RH0*tffcL ++Z 
CS-Q+SWING 
CL-W/QS 

22 CALL TRAJQK2) 

CALL AEROUCALCfNERROR* IGEO.KGPRNT) 

IF (NERR0R.GE.2) RETURN 
CALL TRAJOl(l) 

DRAG«CD*QS 
CALL TRAJ0K2) 

CALL PRQPdCALC.NERROR, IGEO.KGPRNT) 

IF ( NERROR.GE . 2 ) RETURN 
CALL TRAJOl(l) 

TN»THRUST*EN 

TNU-THRSTU*EN 

IF (ABS (QRAG/TN-1.0) .LE.0.03) GO TO 30 
IF (KERR0R.EQ.2) WRITE (6.602) IPHASE. LEG. DRAG. TN 
6C2 FQRHATt/ 17H IN CRUISE. PHASE. 12. 5H. LEG, 12. 46H, THRUST DOES NOT 

1 NATCH ORAL WITHIN 3 PERCENT. /7H DRAG E12. 5, 4X, 4HTN -.E12.5/J 
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15 


2G 


25 


30 


35 


40 


45 


.50 


L55 


L60 


165 


17C 


30 0ELKF-.8444*W*SFC+CD*DELX*SQRT{ 2. *RHO*SWING/ (CL*WN/COSA 

C FUEL USED DURING CRUISE 

WF-WF+DELWF 

DELT-6080.+DELX/VEL 

C CRUISING TIHE TO TARGET 

TIME-TIME+DELT/60. 

W-W-DEL WF 

C TEST FOR TANKS DROP 

WFUSED-WGTQ— W 

IF (IWS.tQ.O) WFUSED-WFUSED-WHISS 
IF (IBS.EQ.O) WFUSED-WFUSED-WBOMBS 
IF (IAS.EQ.C) WFUSED-WFUSED-WAMMUN 
IF (ITS.EO.l.AND.wFUSED.GE.WFEXT) W-W-WETANK 
IF (ITS.EQ.l.AND.WFUSED.GE.WFEXT) ITS-0 
IF (IPRT(IPHASE).EQ.O) GO TO 40 

WRITE (6,603) LEG,RHQ1,RH02,SH,HN,A, VEL, RHO, U, QS> CD, DRAG, TN»SFC, 


603 


40 

604 

50 


C' 



C 

70 


80 

605 

90 


1 DELWF,WF, DEL T.' TIME, W 
FORMAT ( /4X,5KLEG =-I2/3X, 

1 E14 • 7, 5X, 4HHN »E14.7,6X 

2 E14.7,6X,3H0 «E14.7,5X, 

3 E14.7,5X,4HTN -E14.7,4X 

4 E14.7/3X, 6HDELT -E14.7, 
IF (W.GT.WPL) GO TO 50 
IF (KERR0R.EQ.2) WRITE ( 
FORMAT ( //Z7H W.LE.WPL IN 
RETURN 

CONTINUE 

BLTIHE-BLTIMt+TIME/60. 

BLRANG-BLRANG+XCRUSE 

DXCRUS-0.0 

—TEST FOR WEAPONS DROP 
IF ( IW( IPHASE ).EQ«b) GO 
W«W-WMISS 

iws-O 

—TEST FOR BOMBS DROP 
IF (IBIIPHASE) .EQ.O) GO 
W-W-WBOMBS 


6HRH01 -E14.7, 
j3HA -E14.7/4X 
4H0S -F14.7,5X 
, 5HSFC -E3 4. 7, 
3X, 6HTI ME -EI4 

6,604) IPHASE, 
CRUISE. PHASE 


3X,6HRH02 -E14 . 7,5X , 4HSH ■ 
, 5HVEL -E14.7,4X,5HRHD - 
, 4HC0 -E14. 7/3X, 6HDRAG - 
2X, 7HDELWF *E14.7,5X,4HkF 
. 7, 6X, 3HW -E14.7/) 


LEG 

■, 12, 7H 


LEG*, 12/) 


TO 60 


TO 70 


ORIGINAL PAGE IS 
OF POOR QUALITY 


IBS-0 

— TEST FOR AMMO DROP 

IF (IACIPHASE) .EQ.O) GO TO 80 

W-W-WAMMUN 

IAS-0 

IF (IPRINT.EQ.O) GOJO 90 

WRITE ( 6,605) ~W 

F0RHAT(/6X,3HW -E14.7/) 

SMNTl IPHASE) -SHN 
HNT(IPHASE)-HN 
CLT( IPHASE ) -CL 
ALPHAT ( IPHASE ) “ALPHA 
WFT ( IPHASE ) -WF 
TIMETf IPHASE) -TIME 
VELT (IPHASE )-VEL 
SFCT ( IPHASE )-SFC 
TNT( IPHASEJ-TN 
CDT( IPHASE ) -CD 
RLDT( IPHASE )«RLD~ 

WEEL T t IPHASE) *H 


1 

1 
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SUBROUTINE 


175 


18C 


IX);, - T -• j , : : ; i 


CRUISE 76/76 QPT-2 


FTN 4.5+410 


it ART ( I PHASE ) * HTOT 
QT (1PHASE) *0 
SFCUTMPHASE ) “SFCU 
TNUT( IPHASE )»TNU 
CDINSTI IPHASE )*CDINSP 
PRTOTTt IPHASE 1*PRT0T 
XT { I PHASE) *XC RUSE 
RETURN 
END 
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SUBROUTINE ACCEL 


76/76 OPT-2 


FTN 4 . 5+410 


0 


1 

5 


1C 


15 


2C 


25 


30 


35 


4C 


45 


5C 


55 


SUBROUTINE ACCEL ( ICALC# NERROR# IGEQ# KGPRNT) 

REAL MS TART # MEND 

COMMON /TRAJCM/ ALPHA# ARW# CD#CDL# CDO# CL/ DESLF# DRAG#E N#HN# KP» 

1 RLD#SMN#QMAX# RANGE# SFC#SWING#THRUST# TW# ULTLF#W# WAMMUN# WE TANK# 

2 WFUEL# WGTOWT# WHISS# WTOT# HPLWT# MSTART ( 12)#HSTART (12)»TIM(12)# 

3 PSlGT(12)#TD0TSTC12)#NZST(12)#PSIT(12),TD0TITa2)#NZiniZ)» 

4 WBOMBS# WFEXT# SFCU#THRSTU# COINS P# PRTOT #HKFUEL# CRM ACH#FLTO#X( 12)# 

5 FLLAND# TENDUR# BLRANG#BLTINE#TI HT01# TI MTUZ# WFTO# PNAME1 ( 12 ) # 

6 PNAME2 (12 ) #SHNT(12)#HNT(12 )» WFT tl2)#TIHET(12)#XT(12)#CLT(12)# 

7 CDT(12)#ALPHAT(12)#CLIT(12)#CDIT(12)#ALPHITtl2}#WFUSED#WLAND# 

0 QT(12)»RLDT( 12 ) # SFCT< 12> » TNT( 12 ) #CET ( 12) #D Y ( 18 ) # 

9 IAO#IPS#ITS# I WS#IBS# NPHAS E# IDY { 14) 

COMMON /TRAJEX/ ALP1GTU2 J # CDIN.5T ( 12 J # CDiGTt 12 ) # CL1GT ( 12 ) # DECEL# 

1 DEL HP# DEL HP# DXCRUS# FLFAC# FRFURE# HENO( 12)# HENDUR# HMAXP#HMINP# 

2 hPATHT(50I#MEND(1Z)#HENDUR#NT.C12J#PRT0TTC12.).#RADIT(1Z)# 

3 RADST(12)#SF,SFCUT(12>#SMMAXP#SHMINP# 

4 TNUT(12)#VlNDtl2)#VELT(12)#VPATHT( 50) # WART (12) > WC0HBP#WEELT(12)# 

5 WFR ES# WFTOT# WFTOi# WFT02# WFTRAP# WKLAND#XDE5C# XGRLAN#FWGMAX»TOL# 

6 WFIXED# HGCALC# WPL»KGTO# DUMY ( 62 } # I At 12) # IAS# I B( 12) »IBREG# IENDUR# 

7 IMISS#IP(12)#IPHASE#IPLQT»IPRINT#IPRT(12)#IPSIZE#IPST0l#IPST02# 

B IW(12)»IX(lZ)»KERRQR»HILC0H#NCRUSE#NLEGCL#NLEGCR#NLEGLO#NHISS? 

9 NPATH#IDUMY(6) 

DIMENSION ANS14) 

IF (IPRT(IPHASE).NE.O) WRITE (6#6Q0) 

600 FDRHATt ///25H _ ACCELERATION PHASE//) 

IF ( KST ART ( I PHASE j .LT.O.O) SHN-SHNT l IPHASE-1 ) 

IF (MSTART(IPHASE).GE.O.O) SMN»MSTART( IPHASE ) 

IF {HSTARTCIPHASEJ.LT. 0.0) HN«HNT(XPHASE-1> 

IF (HSTART (IPHASE). GE. 0.0) HN-HSTART(IPHASE) 

CALL AT6Z(HN#ANS) 

A»ANS( 4) . ... 

RHO-ANS(l) 

TEMP“1» 8*ANS (3 J 
5QTEMP*SQRT ( TEMP ) 

IPS* IP ( IPHASE ) 

I A0*8 

C COMPUTE INTEGRALS SUHWF AND SUHDT BY TRAPEZOIDAL RULE 

SUHWF-0.0 

SUMDT-0.0 

XACC EL*O.D 

TIMECR*0.0 

PUELCR*0.0 

WF*O.Q 

TiME-0.0 

LEG«0 

IF (IPRT(IPHASE).EQ.O) GO TO 10 

WRITE ( 6# 601 ) LEG#HN#A#RHO#TEMP#SUMWF»SUMDT#SMN#W#IAOf IPS 

601 FORMAT! MX# 5HLEG -I2/5X#4HHN -E14.7# 6X#3HA *£14.7#4X# 

1 5HRH0 *E14.7»3X#6HTEMP -E14.7# 2X# 7HSUMWF -E14.7/2X# 7HSUMDT »# 

2 E14.7# AX# 5HSMN *E14.7#6X»3HW *E14.7#4X#5HIA0 *12# 16X# SHIPS -12/) 
10 IF ( MST ART( IPHASE ) .GE. MEND ( IPHASE ) ) GO TO 120 

IF ( SMNTl IPHASE— 1) . GE . MEND ( IPHASE ) ) GO TO 120 
C — — -START OF LCDP. 

2C LEG-LEG+1 

IF (LEG.GT.50) GO TO BO 
CALL TRAJ0H2) 

CALL PROP ( 1C ALC#NERROR# IGEO#KGPRNT) 


OF POot L nf AGE 13 
^ QUAjjJxy 
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76/76 






\.jv? 


SUBROUTINE ACCEL 


OPT-2 


FTN 4. 5 + 410 05 


60 


65 


70 


75 


BC 


i85 


90 


95 


100 


105 


110 


IF £ NERRUR. GE .2 ) RETURN 
CALL TRAjOKL) 

TN-THRUST*£N 

TNU*THRSTU*EN 

VEL«SMN*A 

30 CL«2.*W/£RH0*SWING+VEL**2) 

CALL TRAJ0K2) 

CALL AERD( ICALC* NERRDR* IGEO*KGPRNT) 

IF (NERRDR. GE, 2) RETURN 
CALL TRAJOl(l) 

DENQH-TN*CDS£ALPHA/57.3)-.5*CD*RH0*SWING*VEL**2 

YWF«SFC*TN/DENOH 

YDT*1 . /DENOH 

IF £ (YWF.GT.O.O) .AND, (YDT.GT.O.O) ) GO TO 50 
DTIME-60. 

DFUEL«TN*SFC*DTIHE/3600. 

W-W-DFUEL 

IF (W.GT.WPL) GO TO 40 

IF tKERR0R.EQ.2J WRITE (6*602) IPHASE* LEG 

602 FORMAT t //26H W.LE.WPL IN ACCEL. PHASE- * 12* 7H LEG»*I2/> 

RETURN 

AO DXCR»VEL*DTlME/6080. 

TIHECR»(7IHECR+DTIME) /60. 

FUELCR-FUELCR+DFUEL 
XACCEL-XACCEL+DXCR 
GO TO 30 

50 IF (LEG.EQ.l) GO TO 60 

SUHWF-SUMwF+.5*(YWF+YWFP)*tSMN-SHNP) 

SUMDT«SUMDT+.5*£YDT+YDTP)*£SMN-SHNP) 

DELX-. 75E5*( YDT+YDTP) *£ SHN— SHNP )*VEL*W*SQTEHP/6080» 
XACCEL-XACCEL+DELX 

WF«(EXP(1.525*SUMWP*SQTEHP/3600.)'-l.)*W 

DELWF-WF-WFP 

W-W-DELWF 

TIME»1.525*SUMDT+SQTEMP*W/60. 

DELT-TIHE-TIMEP 

60 IF ( IPRTtIPHASE).EQ.O) GO TO 70 

WRITE £6,603) LEG* TN* SFC* VEL* SHN* ALPHA * CL* CO* DENOH* YWF* YDT*SUHWF* 

1 DELWF* SUHDT* DELT* DELX* X ACCEL 

603 FOR MAT ( /AX* 5HLEG «I2/5X*4HTN -E1A.7* AX#5HSFC *E14.7*4X* 

1 5HVEL ■E14.7*4X*5HSHN «E14.7*2X*7HALPHA «E14.?/5X,4HCL «E14.7*5X* 

2 AHCD »EI4.7*2X*7H DENOH -E14. 7* 4X* 5HYWP -E1A.7*AX*5HYDT -E14.7/ 

3 2X, 7HSUHWF »E14 .7* ZX*7HDELWF «E1A,7*2X*7HSUHDT -E14.7* 

A 3X* 6H0ELT "E14. 7* 3X* 6HDELX -E14.7/1X* 8HXACCEL -E14.7/) 

C HAKE TEST Ttl SEE IF YOU ARE DONE. IF SO* JUMP OUT DF LOOP. 

70 IF CSHN • GE • HEND 1 1 PHASE) ) GO TO 80 
YWFP-YWF 
YDTP-YOT 
SFfir-SHN 
WFP-WF 
TIHEP-TIHE 

IF £SHN*LEt.S7.QR.SHN»GE.l.l) DELSMN-.05 
IF (SHN.GT.. 87. AND.SHN.LT. 1.1) DELSHN-.02 
SMN-SMN+DELSHN 

IF £ SHN. GE. HEND ( IPHASE ) ) SMN-HENDI IPHASE ) 

C END OF LOOP. GO BACK TO START OF LOOP. 

GO TO 20 
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FTN 4.5+410 


115 


120 


125 


130 


135 


14C 


145 


150 


155 


160 


C FUEL USED DURING ACCELERATION RUN 

80 WF-WF+FUELCR 

C‘ TIME ELAPSED DURING ACCELERATION RUN 

TIME-TI KE+TIMECR 
BLTIHE-BLTIME+TIME/60. 

BLRANG-BLRANG+XACCEL 
IF (W.GT.WPL) GO TO 90 . 

IF ( KERRDR . EQ .2 ) WRITE (6*604) IPHASE 

604 FORMAT { /26H W.LE.WPL IN ACCEL. PHASE-* 

RETURN 

90 VEL«SMN*A 

0«.5*RH0*VEL**2 

DXCRUS-DXCRUS+XACCEL?FI,5AT.UX(.JPfclASEJ) 

C TEST FOR WEAPONS DROP 

IF (IW(IPHASE),EQ.Q> _G0..Ta_lQ.p 

W-W-WMISS 

IWS-0 _ __ 

C — - — TEST FOR BOMBS DROP 
100 IF ( I B ( IPHASE ). EQ.0Jl_GJ3_TJ3_110 
W-W-WB0M8S 
IBS-0 _ 

C TFST FOR AMMO DROP 

110 IF (IA{IPHASE)«EQ.C) GO TO 120 
h-W-WAMMUN 

IAS-0 . ... . 

120 IF (IPRT(IPHASE).EQ.O) GO TO 130 

WRITE (6*6053 SUMWF* SUMDT* TIHE* VEL* a* XACCEL* DXCRUS* WF* W 

605 FORHAT( /2X* 7HSUNWF -E14. 7*2X* 7HSUMDT -E14.7* 3X* 6HTIME - 

1 E14 .7* 4X* 5HVEL «E14.7*6X*3HQ -E14. 7/lX*8HXACCEL -E14.7*1X* 

2 BHDXCRUS -E14. 7*5X*4HWF -E14 .7* 6X*3HW -E14.7/3 

130 SMNTCIPHASEJ-SHN 

HNT (IPHASE 3 -HN 

CLT( IPHASE) -CL 

ALP HAT (IPHASE 3 -ALPHA 

WFT( IPHASE 3 -WF 

TIME T( I PHASE 3 -TIME 

VELTt IPHASE 3 -VEL 

SFCT(IPHASEJ«SFC 

TNT( IPHASE J-TN 

CDT ( IPHASE ) -CD 

RLDT( IPHASE 3 -RLD 

WEELT ( I PHASE 3 -W 

WART(IPHASE3*WT0T 

QT(IPHASE) -0 

SFCUT( IPHASE 3 -SFCU 

TNUT(IPHASE)-TNU 

CDINST (IPHASE 3-CDI NS P 

PRT0TT( IPHASE 3 -PRTOT 

XT(IPHASfc) -XACCEL 

RETURN 

END 
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20 
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40 
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SUBROUTINE COMBAKICALC/NERROR/ IGEOz KGPRNT) 

REAL MSTART/ NT/ NTURNSz LIFT zNZS/NZI/NZST/NZIT 

COMMON /TRAJCM/ ALPHA/ARW/CD/CDL/CDO/CLz DESLF/ DR AG/ ENz HN/ KP/ " 

1 RLD/SMNz 3MAX/ RANG E/SFCz SWING/ THRUST/ TW/ULTLP/W/WAMMUN/ WE TANK? 

2 WFUEL/WGTOWT/WMISS/WTOT/WPLWT/ MSTARTI12) / HST ARK 12) z TIM( 12) / 

3 PS1GK12)/TD0TSK12)/NZSK12)»PSIT(12)/TDQTIK12»/NZIK12)/ 

4 WBOHBS/WFEXT/SFCU/THRSTU/CDINSP/PRTOT/WKFUEL/CRMACH/FLTO/Xt 12)'/ 

5 FLLANO/TENDUR/BLRANG/BLTIME/TIMTO 7IHT02/WFT0/PNAHEK12)/ 

£ PNAME2(12l/SMNU12)/HNTt 12J/WFK 12 j / TIMEK 12 ) /XT { 12 ) /CLT ( 12 ) / 

7 CDT(12)/ALPHAT(12)/CLlTtl2)/CDIT(12)/ALPHIT(12)/WFUSED/WLAND/ 

8 QT(12)/RLDT(12)/SFCTtl2)/TNT(12)/CET(12)/0YClB), 

9 1A0/ IPS/ITS/IWS/IBS/NPHASE/IDY < 14 ) 

COMMON /TRAJEX/ ALP1GT ( 12 ) / CD INST 1 1 2 i / CD1GK 12 1 /CL 1GT 1 12) /DEC EL/ 

1 DtLHP/DELMP/DXCRUS/FLFAC/FRFURE,HEND(12l/HfcNDUR/HHAXP/HMINP/ 

2 HPATHH50)/MENDll2)/MENDURzNK12)/PRTQTK12)/RADIK12)/ 

3 RADST(12)/SF/SFCUT U2I/SMMAXP/SMMINP, 

4 TNUT(12)zVlNDU2)/VELT(12)/VPATHTl50)/WAPK12)zWCOMBP/WEELT(12)/ 

5 WFRES/WFT0T/WFT01/WFTO2/WFTRAP/WKLAN0/XDESC/XGRLAN/FWGHAX/T0L/ 

6 WFIXE0/WGCALC/WPL/WGT0/DUMY{82)/IA(12)/IAS/ IB l 12 ) / I BREG/ IENDUR/ 

7 IMISS/IPC12)/ IfHASE/ I PLOT / I PRI NT/ I PRT ( 12 ) / 1 PS IZE / IPST01/ I PS T02/ 

0 IW(12)/1X(12)/KERROR/MILCOM/NCRUS£/ NLEGCL/NLEGCR/ NLEGLD/ NMI SS/ 

9 NPATH/ IDUHYC6) 

DIMENSION ANS(4) 

HN-HSTARK1PHASE) 

IF (HSTART(IPHASE) .LT.O.O) HN*HNK IPHASE-1) 

CALL AT62IHN/ANS) 

RHO-ANStl) 

A*ANS ( 4 ) 

SMN«M$TART( IPHASE ) 

IF ( MSTART t IPHASE ) . LT.O. 0 ) SMN-S MNTtlP HASE-1 ) 

VEL-SHN+A 

Q«.5*RH0*VEL**2 

OS*Q+SWING 

XCOHB*0.0 

IF (IPRTIIPHASE).NE.O) WRITE (6/600) W/HN/RHOz A/ SMN/ VEL/ Q/ QS / XCOMB 
600 FORMAT! // /19H COMBAT PHASE/// 

1 6X/3HK -E14.7/5X/4HHN -E14.7/4X/5HRH0 *E14» 7/ 6Xz 3HA *E14.7/ 

2 4X/5HSNN *Ei.4i 7/ 4Xz 5HV EL -E14.7/ 6Xz 3H0 -E14.7/ 5X/4HQS -E14.7/ 

3 2X/7HXC0MB -E14.7/1 

C — "•* 

C ONE G FLIGHT 

CL* W /OS __ 

IAQ*B 

CALL TRAJ0K2) 

CALL AEROdCALC/NERROR/IGEO / KGPRNT) 

IF (NERROR.Gt.2) RETURN 
CALL TRAJ0H1) 

DRAG»CD*QS 
IPS* IP ( IPHASE ) 

CALL TRAJ0K2) 

CALL PROPtICALC/NERROR/IGEO/KGPRNT) 

IF (NERR0R.GE.2) RETURN 
CALL TRAJ0K1) 

TN-THRUST+EN _ 

TNU-THRSTU*EN 

PS1G*CTN-DRAG)*VEL/W 
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FTN 4.5+410 0! 


60 


65 


7G 


75 


8C 


85 


90 


95 


LOO 


105 


110 


PS1GT(IPHASE!*PS1G 

ALP1GT { I PHASE ) “ALPHA 

CL1GT(IPHASE)*CL 

COIGKIPHASE J-CD 

IF (IPRT(IPHASE).EQ.O) GO TO 10 

WRITE (6/601) I AO/ CL/ C0» ALPHA/DRAG/ I PS/TN/SPC/PS1G 

601 FORMAT ( /1X/12H0NE G FLIGHT// _4X/ SHI AO -I2/17X/4HCL -E14.7/ 

1 5X/4HCD *E14.7/2X/ 7H ALPHA *E 14. 7/ 3X/ 6HDRAG *£14.7/ 

2 4X/5HIPS * 12/ 17X/ 4HTN “E14.7/4X/5HSFC *E14. 7/ 3X/ 6HPS1G *E14.7/1 

C SUSTAINED TURNS 
C — 

10 NZS-l.OOL 

IFtPSIG.LE.O. ) GO TO 40 

IAO-4 

CALL TRAJ0U2) 

CALL AEROt XCALC/NERRQR/IGEO/KGPRNT) 

IF (NERR0R.GE.2) RETURN 

CALL TRAJOl(l) 

L1FT-CL+QS 

C LOAD FACTOR 

NZS* <LIFT+TN*SIN{ ALPHA/57.2958 > )/w~ 

IF (IPRT(IPHASE) .EQ.O) GO TO 20 
WRITE (6/602) I AO/ CL/CD/ ALPH A/ L IFT/ NZS 

602 FORMATl/lX/lSHSUSTAINED TURNS// 4X/5HIA0 -I2/17X/4HCL -E14,7/_ 

1 5X/4HCD *E14.7/2X/7H ALPHA “E14.7/3X/6HLIFT *E14.7/ 

2 4X/5HNZS *E14. 7/ ) 

20 IF (NZS.LE.DLSLFS GO TO 30 
IAO* 6 

CALL TRAJ01 ( 2 ) 

CALL AERD( IC ALC/ (.ERROR/ IGEO/_KGPRNT) 

IF (NERRQR.GE.2) RETURN 
CALL TRAJOKl) 

NZS-DESLF 

IF (IPRT(IPHASE) .EQ.O) GO TD 30 

WRITE (6/603) IAO/CL/CD/ALPHA/DESLF/NZS 

603 K1RMATC/4X/5HIA0 -I2/17X/4HCL *E14. 7/ 5 X/ 4HCD *E14.7/ 

1 2X/7HALPHA *E14.7/2X/7HDESLF * E14 . 7/4X/ 5HNZS *E14.7/) 

30 IF (NZS. GT. 1.0) GO TO 50 

4C IF (KERR0R.EQ.2) WRITE (6/604) IPHASE/NZS 

604 FORHAT ( / / 36H INSUFFICIENT LIFT IN COMBAT. PHASE*/ I2» 5X/ 

1 4HNZS*/1PE11.4) 

NZS*I.001 

WRITE (6/605! NZS 

605 F0RHAT( /4X/5HNZS -E14.7/) 

50 PHI-ARC0Stl.0/NZ5) 

C TURN RADIUS 

RADlUS*VEL**2/{ 32. 2*TANt PHI) ) 

C TURNING TIME OR NUMBER OF TURNS 
TIHE*TIH( IPHASE ) 

NTURNS-NTCIPHASE) 

IF (TIME. LE. 0.0) TIME-NTURNS*RA0IUS/(9.55*VEL) 

IF (NTURNS.LE.O.O) NTURNS«9.55*TIME*VEL/RADIUS 
XCOHB-VEL*TIME*6O./6O0O. 

BLTIME-BLTIME+TIME/60. 

BLRANG*BLRANG+XCOMB 
C FULL USED DURING TURNS 
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0 


115 


LHC 


L25 


L3C 


135 


140 


145 


150 


155 


16C 


165 


17C 


WF-SFC*TN*TIME/60. 

C COMBAT ENERGY 

CE»3600.*WF*PS1G/(SFC*TN) 

TDQTS*57i3*V£L/RADIUS 

RADS-RADIUS 

IF (iPRTtlPHASEKEQ.O) GO TO 60 

WRITE (6,606) PHI, RADIUS, TIME, NTURNS, W F, TDGTS, XC OMB 
6C6 FQRNAT(/4X,5HPHI » E14 . 7, IX, 8HRADIUS *E14 .7, 3 X, 6HT IMF *E14.7» 

1 IX, BHNTURNS «E14.7,5X,4HWF -E14.7/2X, 7HTD0TS *E14.7, 

2 2X, 7HXC0MB -E14.7/) 

6C SMNT( IPHASE )«SHN 

HNT ( IPH AS6 ) «HN 
CLT(IPHASE) "CL 
ALPHATUPHASE)-ALPHA 
WFTl IPHASE ) *WF 
VELT(IPHASE)*VEL 
TIHETI I PHASE) “TIME 
SFCT ( IPHAS E ) «SFC 
TNT( IPHASE )*TN _ 

CDT (IPHASE ) *CD 
RLDTdPHASE )*RLO 
WART (IPHASE)»WTQ7 
OT ( IPHASE ) “Q 
SFCUT(IPHASE)«SFCU 
TNUT (IPHASE )*TNU 
CDINST(IPHASE)«COINSP 
PRTOTTCIPHA$El«PRTOT 
XT(IPHASE)*XCOMB 
TDOTST( IPHASE )*TDOTS 
RADST ( IPHASE ) *RADS 
NZST(IPHASE)«NZS 
CET (IPHASE )*CE 

C 

C MAXIMUM INSTANTANEOUS CONDITIONS 

IAQ-5 ' 

CALL TRAJQK2) 

CALL AERO ( ICALC, NERROR, IGEO, KG PR NT 1 
IF (NERROR. GE. 2) RETURN 
CALL TRAJOl(l) 

CLHAX»CL 

LIFT«CLMAX*QS 
C LOAD FACTOR 

NZIM LIFT+TN+SIN ( ALPHA/57. 2958 1 ) /W 
IF (IPRT( IPHASE) .EQ.O) GO TO 70 
WRITE (6,607) I AO, CL, CD, ALPHA, L 1FT , NZI 
607 FORMAT ( /IX, 19H INSTANTANEOUS TURNS// 4X,5HIA0 *12, 17X,4HCL -E14.7, 

1 5X,4HC0 «E1'-,7,2X,7HAIPHA »E14.7,3X,6HLIFT *£14.7/ 

2 AX, 5HNZI -E14.7/) 

70 IF (NZI.LE.DESLF) GO TO 80 
IA0*6 

CALL TRAJ0K2) 

CALL AERO( ICALC, NERROR, I GEO, KG PR NT) 

IF (NERROR. GE. 2) RETURN 
CALL TRAJ0H1) 

NZI-DESLF 

IF (IPRTUPHASEl.EQ.O) GO TQ 80 
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20C 


205 


210 


215 
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WRITE ( 6, 608 ) IAFI,CL,CD, ALPHA,DESLF, NZI 
608 FORMATC / 4X,5HIAQ -IZ,17X,4HCL -E14. 7, 5X, 4HCD -E14.7, 

I 2X,7HALPHA «E14.7,2X,7HDESLF -E14. 7/4X, 5HNZI -E14.7/) 

80 IF CNZI.GT.1.01 GO TO 90 

IF <KERROR.EQ.25 WRITE (6,609) IPHASE, NZI 
6C9 FGRHAT ( / / 36H INSUFFICIENT LIFT IN COMBAT. PHASE-, 12, 5X, 

1 4HNZI-, IP Ell .45 
RETURN 

C PS FDR MAXIMUM INSTANTANEOUS TURN RATE 

90 DRAG«CD*QS ORIGINAL PAGE IS 

PSI«£TN*C0$(AlPHA/57.3)-DRAG)*VEL/W OF POOR DTrArmrS 

ALPHI-ALPHA W&hLl X 

CLI-CL. _ 

COI-CD 

PHI-ARC0S11.0/NZI) 

C MAXIMUM INSTANTANEOUS TURN RATE ANO RADIUS 
RADIUS- VEL**2/ <32. 2*T AN (PHI I ? 

TDOTI-57. 3* VEL /RADIUS 

RADI-RADJUS 

ynK-yp 

IF (W.GT.WPL) GO TO 100 
IF (KERROR.EO.2J WRITE (6,610) IPHASE 

610 F0RHAT(//27H W.LE.WPL IN COMBAT . PHASE-,12) 

RETURN 

C TEST FOR WEAPONS DROP 
ICO IF (IW(IPHASE).EQ.O) GO TO 110 
W-W-WMISS 
IWS-0 

C TEST FOR BOMBS DROP 
110 IF tlB(IPHASE).EQ.O) GO TD 120 
W'-W-WBOMBS 
IBS-0 

C TEST FOR AMMO DROP 

12C IF (IA(IPHASb).EQ.O) GO TO 130 

W-W-WAMHUN 
IAS-0 

130 IF (IPRT(IPHASE).EQ.O) GO TO 140 

WRITE (6,611) DRAG, PSI, PHI, RADIUS, TDOTI,W 

611 FCRHAT(/3X, fcHDRAG -E14. 7,4X, 5HPSI « E14‘.7,4X, 5HPHI -E14.7, 

1 IX, 8HRADIU5 -E14. 7, 2X, 7HTD0TI -E14.7//6X,3HW -E14.7/) 

140 WEELTUPHASFl-W 
PSITUPHASEJ-PSI 
TDOTITt IPHASE )-TDOTI 
RADITCIPHASE)-RADI 
NZIT(IPHASE)-NZI 
ALPHITt IPHASE J-ALPHI 
CL IT ( IPHASE )-CL I 
CDIT(IPHASE)-CDI 
RETURN 
END 
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SUBROUTINE LOITER{ICALC»NERROR>IG£G>KGPRnT) 

REAL MSTART>KP>NZLUIT 

COMMON /TRAJCM/ ALPHA>ARW>CD>CDL>CDa>CL>DESLF>DRAG>EN>HN> kp> 

1 R LD>SMN>GMAX>RANGfc*SFC> SWING* THRUST# TW>ULTLF> W> WAHMUN> WE TANK, 

2 wFUEL>W£T0WT>WMISS*tfTDT»WPLUT*MSTART{12).>HSTARTC12)*TIM( 12) > 

3 PSlGT(12i»TD0TST(12}>NZST{12)>PSIT(12)>TD0TIT£12)>NZmi2]> 

A W8DM8S*WF£XT>SFCU*THRSTU*CQINSP>PRT0T>WKFUEL>CRHACHjiFLTQ>X(12)> 

5 FLLAND>T£n0UR> BLRANG* BLTINE* TIMT J1> TIMT02»WFTG»PNAME1(12)> 

6 PNAhE2£1Z}>SHNT(12>>HNT(12)>WFTU2J>TI!1ET(12)>XT£12 )>CLT£12)> 

7 CDT(12)>Al PHATf I2}»CLIT(12)»C0IT(12)*ALPHIT£12)>WFUSED*WLAND* 

8 CT(12) ,RLDT(12>,SFCT(12)>TNT{12)>C£T{12)>DY(18}, 

9 IA0>IPS>ITS>IWS>I BS> NPHAS£> IDY ( 14) 

CUMMUN / TRAJEX/ AL P1GTU2) >CDINS T ( 12 ) * CD1GT( 12 ) > CL 1GT (12 ) > DECE L* 

1 D EL HP* DELMP* DXCRUS* FLF AC>FRFUR £> HEND ( 12 ) > HE NOOR* HMAX P>HM I NP * 

2 HPATHTI50)*MbND£12)»MENDUR*NT£12)>PRTQTTll2)>RADIT£12)> 

3 RADST(12)>SP*SFCUTC12) * SMMAXP* SMHINP> 

A TNUT£12)>VIND£12),V/£LT(12)>VPATHT£ 50)>WART<12)> WC0MBP> WE £ LT ( 12 ) , 

5 WFRES>rfFT0T>WFTOl>WFTO2*WFTRAP>WKLAND>XD£SC>XGRLAN>FWGHAX>T0L> 

6 WFIXED>WGCALC>WPL>WGTQ> OUMY (82)>IA(12)>1AS> IB (12 ) > I BR EG* I ENDUR* 

7 IHISS*IP£12)*IPHASE*IPLOT»IPRINT*IPRT£12)*IPS1ZE*IPST01*IPST02* 

8 lH(12)>lX(12)>KERR0R>MILCaM>NCRUSE>NL£GCL>NLEGCR>NLEGL0,NMISS* 

9 NPATH> IDUMY(6) 

DIMENSION A NS (4)*HNTAB(10)*SHTAB ( 10) 

DATA HNTAfl/0.> 5000. >10000. >15000. *20000. >25000. >30000.* 

1 35000, >37000. >100000./ 

DATA SHT AB/ 34500. >32900. >31900. >30500. >29400. >28 200. > 27110. > 

1 25950. *209C0. >20900./ 

DATA DMQN/G. 0/ 

IF UPRTtlPHASEJ.NE.O) WRITE (6>600! 

FORMAT £ ///19H LOITER PHASE//) 

IPS-4 

SMNIN*MSYART( IPHASE ) 

IF £ MSTARTUPHASE ) .LT.O.O) SHNIN-SMNT£ IPHASE-D 

SMN-SMNIN 

HN*HSTART( IPHASE) 

IF (HSTARTIIPHASE) .LT.O.O) HN-HNK I PHASE-1 > 

If CIENOUR.EQ.O.DR.HSTARTCIPHASEl.HE.O.O) GO TO 10 
CALL G0LDEN£200C.> 65000 •> 2> 12*XX»Y>NbRRQR>IGE0>KGPRNT) 

IF (NERR0R.GE.2) RETURN 
HN-XX 

CALL AT 62(HN* ANS ) 

RHO-ANS(l) 

A*ANS £4 ) __ 

VEL-SMN*A 

IF £ SMN.GT .0.0) CL*2.*w/ £ RHO*SW ING+VEL **2 ) 

TIHE*TIMt IPHASE ) 

NLEGS*TIME/fc0.+1.0 
IF (NLEGS.LT.l) NLEGS-1 
IF (NLEGS.GT.6) NLEGS-6 
IF (NLEGLO.GT.O) NLEGS-HLEGLO 
DELT»TIME/FL0AT£NLEGS) 

WP*C.O 
XLOIT *0 .0 

IF tlPRNI PHASE) . EO .0 ) GO TD 20 

WRITE ( 6> 6G1 ) SMNlN*SMN*HN>RHO>A*tfEL»CL*TIME*NLEGS*DELT>WF>XLQIT 
FORMAT! /2X> 7HSMNIN »E 14. 7, 4X> 5HSHN *E14.7>5X>4HHN - 
1 £14 • 7> 4X> 5HRH0 *614* 7* 6X> 3HA -E14, 7/4X> 5HVEL -E14.7* 5X»4HCL -> 
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2 fcl4.7*3X*6HTIHE “E14 .7»2X*7HNLtGS *I2*15X*6HDELT -E14.7/5X* 

3 4HMF »E14.7*2X*7HXLUIT -E14.7/) 

20 DO 90 LEG*1,NLEGS 

1+ (SMN1N. GT« 0*0) GO TO 30 

C* PERFORM OPTIMUM LOITER MACH NO. SEARCH. 

CALL GQLD£U(.3*.9*3»8>XX*Y*N£RR0R* IGEO j KGPRNT) 

IF (NERRGft.GE.2) RETURN 

SMN-XX 

VEL«SMN*A 

IF t LEG.EQ ,1 ) CL0UM«2.*W/(RH0*SWING*VEL+*2> 

CL-CLDUH 
GO TO 40 

C GPTION TO INCREASE ALTITUDE AS A/C WEIGHT DECREASES* MACH KEPT CONST 

30 RH01-RHQ 

RH02*2.*W/(CL*SWING*VEL**2) 

CALL TAINTIHNTAB* SNTAB*HN* SN* 10* 1*NERR *DMQN) 

HN-HN+2 . +SH* ( RH01-RH02) / C RH01+RH02) 

CALL AT62(HN#AN$) 

A* ANS (4 5 

VEt«SMN*A 

RH0«2.*W/(CL*SWING*VEL**2) 

C END OF INCREASE ALTITUDE OPTION 

40 IA0*8 

CALL TRAJ01121 ™ T nnT*T nAnr. to 

CALL AERO( ICALC*NERROR* IGEO* KGPRNT) UKitriMAL FAGL lb 

IF {NERR0R.GE.21 RETURN OF POOR QUALITY 

CALL TRAJ0K1J 
DRAG*.5*RH0*SWING*CD*VEL+*2 
CALL TRAJ01121 

CALL PROP ( ICALC* N ERROR* IGEO* KGPRNT! 

IF (NERR0R.GE.2) RETURN 
CALL TRAJOKIJ 
TN»THRU5T*EN 
TNU-THRSTU+EN 

C CHECK LOITER RADIUS. 

RALDIT»X(IPHASE) 

IF (RALOIT.Lc.O.O) GO TO 60 

NZLOIT* t «5*CL *R HO*SWI NG+VEL *+2+TNi'S I N( ALPHA/ 57.3)1 / W 
IF (NZLOIT.LT. 1.30001) NZL0IT*1 .00001 
C0SPHI«1,/NZL0IT 

C —RADIUS FOR OPTIMUM LOITER. 

R10IT»{5„UE-6*VEL**2*C0SPHI>/SQRTC l.-C0SPHI**2) 

IF ULOIT.LE.RALOIT) GO TO 60 

DO 50 J»l* 3 

DRAGP-DRAG 

CL“ (W/( SWING*. 5*RH0+VEL**2) )*SQRT(l.+( 5. 11E-6*VEL**2/R ALO I T ) **2 ) 
IA0»8 

CALL TRAJ0K2) 

CALL AERO(ICALC*NERROR*IGEO*KGPRNT) 

IF (NERR0R.GE.2) RETURN 
CALL TRAJOl(l) 

DRAG=.5*KHa*$WlNG*CD*VEL**2 

IF ( ABS tDRAG/ DRAGP—l *0 ) .LE.0.03) GO TO 60 

IF (J.EQ.3) WRITE (6*602) IPHASE*LEG 

602 F0RMAT(//55H INSUFFICIENT THRUST FOR REQUIRED LOITER RADIUS. PHASE 
1«*I2*7H LEG** 12/ ) 

CALL TRAJ01C2) 
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160 


165 
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CALL PROPCICALCf NtkROR#IGLQ,KGPRNT) 

If (NERR0R.GE.2) RETURN 
CALL TRAJOl(i) 

TN»TrtRU3T*£N 
TMj-THRSTU+EN 
50 CONTINUE 

60 TEST=DRAG/TN 

If ( (ABS( TE5T-1.0) .Lt. 0.05). AND. (SHNlN.LE.O.U) > GO TO 70 

IF ( TEST.GT .1.0) TN-DRAG 

CD-2.*TN/(t<H()*SWlNG*VtL**2) 

DF AG«.5*RHCj*SwING*CO*VEL**2 

IF (KtRRUR.EO.2) WRITE (6,603) I PHA$ E, LEG, TE 5T 

603 FORMAT (//’7F> IN LOITER, PHASE, 12, 5ri, LEG,I2,8H, DRAG «,Fb.3» 

1 7H OF TN./26H CO CHANGED TO HATCH THRUST./) 

70 DELWF*SFC*TN*0ELT/60. 

WF-wF+OtLwF 
W-W—OE L Wf 

QELX«VEL*DELT*60./6J60. 

XLOIT-XlUIT+OELX 

IF C IPRT( IPHASE ) .£0.0) GO TO 80 

hhlTc (6,604) LEG,SHNIN,SHN,ALPHA,CL,C00,C0,KP,RH01,RH02, SH, HN, A, 

1 FhQ,VtL,0RAG, TN,DELWF,WF,W,OLLX, XL 01 T, SFC ,NZ LOI T, COSPHI, R LOI T 

604 FLiRMT(/4X,SHLfcG -I2/2X, 7HSMNIN *tl4.7,4X, 5HSMN -E14.7, 

1 2X» 7HALPHA * E14 . 7, 5X , 4HCL » E 14 . 7, 4X, 5HCD0 *■£ 14. 7/ 5X , 4HC0 -E14.7, 

2 5X, 4HKP «fcl4.7,iX,6HRH0l «E 14. 7, 3X, 6HRH02 -E14. 7, 5X, 4HSrt -E14.7/ 

3 5X, 4HHN -E14.7, 6X,3HA «E 14 » 7, 4X, 5HRH0 -tl4. 7»4X, 5HVEL -E14.7, 

4 3X,6HDRAG «E14.7/5X,4HTN -E14. 7, 2X, 7H0ELWF - E 1 4. 7, 5X , 4HWF -E14.7, 

5 6X r 3HW «E14.7,3X,6H0ELX »E14.7/2X,7HXLCIT «E14. 7,4X, 5HSFC «E14.7, 

6 IX 8HNZL0IT « E14. 7, 1 X , 8HC0S PH I «£14 .7,2X, 7HRLCI T -E14.7/) 

8C IF (W.GT.wPL) GO TO 90 

IF (KERROR.Eu.2) WRITE (6,605) IPHASE, LEG 
6C5 FORMAT ( / /27H W.LE.WPL IN LOITER. PHASE-, 12, 7H LEG-,I2/» 

RETURN 

90 CONTINUE 

Q-.5*k'H0*VEL**2 
BLTIHE.-8LT1HE + TIME/60. 

BLRANG- 8!. P.ANG + XLO I T 

C TEST FOR WhAPQNS DROP 

IF (Ih(IPMASfc).EQ.O) GO TO 100 

W-W-WMISS 

IUS-0 

C TEST FOR 90MBS DROP 

ICO IF (IB(IPHASE).EO.O) GO TO 110 

W-W-WBOMBS 
IBS-0 

C TEST FOR AMMO OROP 

110 IF C 1 A t IPHASE). EQ.O) GO TO 120 

W-W-WAflMUN 
1AS-0 

120 SUNT (IPHASE )-SMN 
HNT ( IPHASE)-HN 
CLT( IPHASE ) -CL 
ALPHAT( IPHAS£)-ALPHA 
WFT ( IPHASE ) -W F 
TIHET(IPHASE)*TIHE 
VELT(IPHASE)-VEL 
SFCTl IPHASEl-SFC 
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TNT ( IPHASE ) *TN 

CDTt IPHASE ) *CD 

RLDT(IPHASE)»RLO 

WtELT (IPHASE ) “W 

WART (IPHASE J-HTOT 

QT( IPHASE) *0 

SFCUT f I PHAS E ) *S FCU . 

TNUTCIPHASE)«TNU 

CDINST U PHASE J-COINSP 

PRTOTT ( IPHASE ) “PRTQT 

XTCIPHASE)*XLOIT 

RETURN 

END 


ORIGINAL PAGE IS 
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I SUBROUTINE DESCNTl ICALC, NERROR, IGEO,KGPRNT) 

REAL KP 

COMMON /TRAJCM/ ALPHA, ARW, CD# CDL,CDO,CL, DESLF, DRAG, EN, HN, KP, 

1 RLD,SNN,aMAX, RANGE, SFC, SWING, THRUST, TW,ULTLF,W,WAMMUN, WETANK, 

5 2 WFUEL, WGTOWT, WNISS, WTQT, WPLWT, MSTART(12), HSTART (1ZJ,TIM(1Z), 

3 PSlGT( 12 ),TDaTST(lZ),NZ$TCi 2 >,PSITU 2 ),TDaTIT( 12 ),NZIT( 12 ), 

4 WB0MBS,WFEXT,SFCU,THRSTU,CDINSP,PRT0T,WKFUEL,CRMACH,FLTD,X(12), 

5 FLLAND,TENDUR,BLRANG,BLTIME,TIHT01,TIHT02,WFTD,PNAME1U2), 

fc PNAME2 (12 ) ,SMNT (12)»HNT(12),WFT{12),riMET(12),XT(12),CLT{12)» 

1C 7 CDT(12),ALPHAU12),CHT(12),CDIT(12), ALPHIT (12 > , WFUS£D,WL AND, 

8 QT(12),RLDT(12),SFCT(12),TNT(12),CET(12),DY(18}, 

9 IaC,IPS,ITS,IWS,IBS,NPHASE,IDY(14> 

COHMON /TRAJEX/ ALP1GT(12 ), CDINST (12 ), CQ1GF (12) ,CL1GT(12) , DECEL, 

1 DELHP,DELMP, DXCRUS, FLFAC, FRFURE,HEND ( 12) ,HENDUR, HMAXP,HM I NP , 

15 2 HPATHT (50>,MEND(12) ,MENDUR,NT( 12), PRTQTTC12 ) ,RADIT( 12), 

3 PADST(12),SF,SFCUT(12),SMMAXP,SMMINP, 

4 TNUT(12),VIND(12),VELT(12),VPATHT(50) *WART(12),WCOMBP,WEELT(12), 

5 kFRES,WFTOT, WFT01,WFT02, WFTRAP, WKL AND, XDESC, XGRLAN, FWGMAX,TQL, 

6 VFIXED,WGCALC,WPL,WGT0,DUMY(82),IA(12),IAS,IB{12),IBREG, IENDUR, 

2C 7 IMISS,IP(12),IPHASE,IPL0T,IPRINT,IPRT(12),IPSIZE»IPSTQ1, IPST02, 

6 IW(1Z),IX(12)»KERROR,MILCOK,NCRUSE,NLEGCL,NLEGCR,NLEGLO,NHISS, 

9 NPATH,IDUMY(6) 

DIMENSION ANS (4) 

IF (IPRT(IPHASE).NE.O) WRITE (6,600) 

25 600 FDRMAT(///20H DESCENT PHASE//) 

HNINT«HNT(IPHASE-*1) 

HNFIN*HENO( IPHASE) 

IPS-5 

VEL* VEl T( IPHASE— 1 ) 

3C NLEGS-5 

DbLH*( HNFIN—HNINT ) /FLOAT (NLEGS ) 

HN=HN2NT 
TIME *0*0 
WF«0.0 

35 XDESC-0.0 

IF (IPRT(IPHASE).EG.O) GO TO 10 

WRITE { 6, 6C1 ) HNINT,HNFIN,CL,RLD,VEL,HN,DELT, TIME, XDESC, WF 
6C1 F ORMAT C /2X» 7HHNINT =E14.7,2X,7HHNFIN « E14. 7, 5X, 4HCL - 

1 E14.7,4X,5HRLD *E14.7,4X,5HVEL »E14.7/5X, 4HHN “E14. 7, 3X» 6HDELT 
4C 2 E14»7,3X,6HTIME -E14 .7, ZX, 7HXDESC *E14.7,5X,4HWF «E14.7/ ) 

10 DO 30 LEG*1, NLEGS 

CALL AT62 (HN, ANS ) 

KHO*ANS ( 1) 

A*ANS(4 ) 

45 SMN“VEL/A 

CALL TRAJ3K2) 

CALL PROPt ICALC, NERROR, IGEO,KGPRNT) 

IF (NERROR, .,2) RETURN 
CALL TRAJOIU) 

5C TMTHRUST+EN 

TNU»THRSTU*EN 
I A0*9 

CALL TRAJ0K2) 

CALL AERO (ICALC, NERROR, IGEO,KGPRNT) 

55 IF (NERROR. GE. 2) RETURN 

CALL TRAJOICI) 

CL*SQRT(CDU/KP > 
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65 
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75 


80 


85 


90 


95 


IOC 


105 


11C 


602 


20 

603 

30 


SING-2. *COO/CL 
C 0 SG-SQRTC 1 .-SINGA +21 
VEL-S«RT( 2 .*W*CQSG/{CL+SWING*RH 0 ) ) 

D£LT— DELH/(VEL*SING) 

TIME-TIME+DELT/ 60 . 

DELWF-SFC+TN*DELT/ 3600 . 

WF*WF+DELdF 

DELX-VEL*DELT*C 0 SG/ 608 G. 

XDESC- XQ ESC+DELX 
HN-HN+DELH 
W-W—DELWF 

IF UPRT(IPHASE).EQ.O) GO TO 20 
WRITE C 6 , 602 ) LEG, RHO, A,S MN, TN, SFC, DEL V, VEL, COSG, SING, DELT, TIME, 

1 DELX,XDESC,HN,CL,CD,KP,OELWF,WF,W 
FGRHAT(/AX, 5 HLEG - 12 /AX, 5 HRH 0 - E 1 A, 7 , 6 X, 3 HA -E 1 A. 7 »AX, 

1 5 HSMN -E 1 A. 7 , 5 X, AHTN -E 1 A. 7 , AX, 5 HSFC -ElA. 7 / 3 X, 6 HDELV «E 1 A. 7 ,AX, 

2 5 HVcL -E 1 A. 7 , 3 X, 6 HC 0 SG « E 1 A. 7 , 3 X, 6 HSING -E 1 A. 7 , 3 X, 6 HDELT -E 1 A. 7 / 

3 3 X, 6 HT IME -E 1 A. 7 , 3 X, 6 HDELX -E 1 A. 7 , 2 X, 7 HXDESC -E 1 A. 7 , 5 X,AHHN *, 

A E 1 A. 7 , 5 X, AHCL -E 1 A. 7 / 5 X, AHCD -E 1 A. 7 , 5 X,AHKP -E 1 A. 7 , 2 X, 7 HDELWF 
5 ElA. 7 , 5 X,AHwF *E 1 A. 7 , 6 X, 3 HW -E 1 A. 7 /) 

IF CW.GT.WPL) GO TO 30 

IF (KERR 0 R.EQ. 2 ) WRITE ( 6 , 603 ) IPHASE, LEG 

FORMAT ( // 28 H W.LE.WPL IN DESCENT. PHASE-, I 2 , 7 H LEG-, 12 /) 

RETURN 

CONTINUE 

CALL AT 621 HN, ANSI 
A-ANS(A) 

SHN-VEL/A 


PAGE Is 
of POOB SffAury 


RHO-ANS ( 1 ) 

0-.5*RH0*VEL+*2 

BLTIHt-BLTIKE+TIME/60. 

BLRANG-BLRANG+XDESC 

C TEST FDR WEAPONS DROP 

IF (IWIIPHA5EKEQ.0) GO TO AO 
W*W— WHI SS 

iws-o 

C TEST FOR BOMBS DROP 

AC IF (IB(IPHASE).EO.O) GO TO 50 

W-W-WBQMBS 
IBS-0 

C TEST FOR AMMO DROP 

50 IF (lA(IPHASE).EQ.O) GO TO 60 

W-W-WANMUN 
IAS-0 

60 IF (IPRT(IPHASE).EQ.O) GO TO 70 

WRITE ( 6, 60A ) HN, A, SHN, RHO, Q, DXCRUS, W 
60 A FORHAT( /5X,AHNN -E1A. 7, 6X, 3HA - E1A. 7, AX, 5HSHN -E1A.7, 

1 AX, 5HRH0 *E1A.7,6X,3HQ «E1A.7/1X,0HDXCRUS -E1A.7,6X,3HW -E1A.7/) 
70 SUNT ( IPHASE ) -SUN 

HNTC IPHASE) *HN 
CLT(IPHASE)*CL 
ALPHAT ( IPHASE)- ALPHA 
WFT (IPHASE ) -WF 
TIHET( IPHASE ) -TINE 

VELT(IPHASE)-JfEL 

SFCTdPHASE ) -SFC *" 

TNTUPHASE)-TN ... ^ 
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COT( IPHASE ) *CD 
RLDTUPHASE ) = RLD 
WEELTt IPHASE J *W 
WART ( IPHASE >WTOT 
QT(IPHASEi*Q 
SFCUTt IPHASE )*SFCU 
TNUTl IPHASE )-TNU 
CDINSK IPHASE ) “CDINSP 
PRTDTTt IPHASE )*PRTQT 
XTIIPHASEJ-XOESC 
RETURN 
END 
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SUBROUTINE PA THS ( I CALC,NERR 0 R, IGEO, KGP RNT J 
PEAL KP 

COMMON /TRAJCM/ ALPHA, AR W, CD, CD L, CDO,C L,DESLF,DR AG,EN, HN, KP, 

1 RLD,SMN, QMAX, RANGE, SFC, SWING, THRUST, TW,UlTLF,W,WAMMUN, WE TANK, 

2 WFUEL,WGTOWT, WHIS 5 ,WT 0 T,WPLWT,HSTART( 12 ), HS TART ( 12 ),TIM£ 12 ), 

3 PS 1 GT(I 2 ),TD 0 TST( 12 ),NZST( 12 ),PSIT( 12 ),TD 0 TIT( 12 ),NZIT( 12 ), 

A WB 0 MBS,WFEXT,SFCU,THRSTU,CDINSP,PRTOT,WKFUEL,CRMACH,FLTU,X( 12 ), 

5 FLIAND,TENDUR,BLRANG,BITINE,TIMT 01 ,TIMI 02 ,WFT 0 >PNAME 1 ( 12 ), 

6 PNAME 2 tl 2 ),SMNT( 12 ),HNT{ 12 ),WFTtl 2 f,TIMET{ 12 ),XT( 12 ),CLT{ 12 ), 

7 CDT( 12 ),ALPHATU 2 ),CLITf 12 ) ,CDlTt 12 ) , ALPHITU 2 ) »WFUSED,WLAND, 

8 QT( 12 ),RLDTC 12 ),SFCT( 12 ),TNTtl 2 ),CET( 12 ),DT< 10 ), 

9 1 A 0 ,IPS,ITS,IWS,IBS,NPHASE,IDY( 1 ^) 

COMMON /TRAJEX/ AL PIGTf 12 ) , CDINSTU 2 ) , CDlGTt 12 ) , CL 1 GTC 12 ) , DECEL, 

1 DtLHP,DELMP, 0 XCRUS,FLFAC,FRFURE,HEND( 12 ),HENDUR,HMAXP,HHINP, 

2 HPATHT( 50 ),MENDC 12 ),MENDUR,NT( 12 ),PRTOTT( 12 J,RADITtl 2 ), 

3 PAD 5 T ( 12 ) ,SF,SFCUTtl 2 ),SMMAXP,SMMINP, 

<t TNUT{ 12 ),VINDC 12 ),VELTU 2 ),VPATHTI 50 ),WART( 12 ),WC 0 MBP,WEELT{ 12 ), 
5 WFRES,WFTDT,WFT 01 ,WFT 02 ,WFTRAP, WKL AND, XDESC,XGRL AN, FWGMAX,TOL, 
fc WFIXED,WGCAlC, WPL,WGT 0 ,DUMY(B 2 ),IAU 2 ),IAS, IB (12 ), I BREG, IENOUR, 

7 IMISS,IP( 12 ) , IPHA 5 E,IPL 0 T,IPRINT,IPRT( 12 ),IPSIZE,IPST 01 , IPST 02 , 

0 IW( 12 ), 1 X 112 ), KERRQR,MILC ON, NC RUSE, NL EGCL, NLEGCR, NLEGLQ, NMI SS, 

9 NPATH, IDUMY t 6 ) 

DIMENSION ANS 14 ) 

IPS*IP { IPHASE ) 

T IME* 0 « 0 
XP ATH* 0 « 0 
WF- 0.0 

IF (IPRT(IFHASE) .NE.O) WRITE ( 6 , 600 ) IPS,TIME,XPATH,WF,W 
600 FORMAT ( /// 17 H PATH PHASfc/// 4 X, 5 HIPS *I 2 , 15 X, 6 HTIME -E 1 A. 7 , 

1 2 X, 7 HXPATH » ElA • 7 , 5 X, 4 HWF -El* . 7 , 6 X, 3 HW "E 1 A. 7 /J 


ORIGINAL’ PAGE IS 
OP POOR QUALITY 


C— MAKE TEST TO SEE WHETHER VM IS MACH OR VELOCITY. 

IF CVM.GT.1C.0) GO TO 10 
SMN1-VM1 
SMN2»VM2 
SMN-VM 
VEL1»SMN1*A 
VEL2*SMN2*A 
VEL-SHN*A 
GC TO 20 
1C VEL1-YM1 


VEL2-VM2 1 
VEL-VM 1 
SMNlaVELl/A 1 
SMN2-VEL2/A 

SMN* VEL / A 1 


DC 60 I»2, NPATH 
LEG»I-1 

H1«HPATHT(I-1) 

H2-HPATHTU) 

HN«.5*(H1+H2> 

DELH*H2-H1 

CALL AT62 (HN,ANS ) 

RhO*ANS(l) 

A*ANS (<♦ ) 
VM1*VPATHT( 1-1) 
VM2*VPATHT (I ) 
VM-.5*{ VM1+VM2) 


57 


• I 

: Ir 




i 




SUBROUTINE PATHS 76/76 OPT-2 FTN 4.5 + 410 05 


20 DELV-VEL2-VEL1 

C CALL PROP TO GET TN* SFC. 

60 CALL TRAJ01UJ 

CALL PROP(ICALC*NERROR* IGEO*KGPRNT ) 
IF ( NERROR .GE .2 ) RETURN 
CALL TRAJ01C1) 

TN-THRUST+EN 

65 TNU-THRSTU* EN 

C CALL AERO TO GET CDO* KP. 

IAC1-1 
ALPHA-5.0 
CALL TRAJ0112) 

7C CALL AEKO(ICALC*NERROR»IGEO*KGPRNT) 

IF {NERROR, GE. 2) RETURN 
CALL TRAJOl(l) 

Q=.5*RH0*VEL**2 

QS-O+SWING 

75 Cl-as/w 

C2«VEL*0ELV/C32.2*DELH)+l.O 

C3-TN/W-C1+CD0 


C4-CC2/KP)**2 

C5-SQRT ( .25*C4**2-C4*C3 / £ Cl+KP ) +C4/C 1**2 ) 

8C C CALCULATE FIRST APPROXIMATION TO CL. 

CL1«SQRTI-.5*C4+C3/(C1*KP)+C5) 

C CALL AERO TO GET CD* ALPHA. 

CL-CL1 
I AQ-8 

85 CALL TRAJ01 (2 ) 

CALL AERCJ< ICALC* NERROR* IGEQ* KGPRNT) 

IF (NERROR. GE. 2) RETURN 
CALL TRAJOKI) 

Cfc-( (TN/W)*C0S(ALPHA/57.3)-Cl*CD!/C2 

9C C CALCULATE SECOND APPROXIMATION TO CL. 

CL2-SQRT(1,-C6**2)/C1 
C”"'—”C ALL AERU TD GET CD* ALPHA, 

CL-CL2 

CALL TRAJ01C2) 

95 CALL AERO (ICALC* NERROR* I GEO* KGPRNT) 

IF (NERROR. GE. 2) RETURN 
CALL TRAJOKI) 

C CALCULATE FLIGHT PATH ANGLE. 

$lNG«((TN/W)*COS(ALPHA/57»3)-Cl*CO) /C2 
LOG IF ( (SING. GE, 0.0). OR, (DELH.LE.O.O) ) GO TO 30 

NERRUR-2 

WRITE (6*601) IPHASE* LEG 

601 FORHATt /////33H *****FATAl ERROR IN PATH. PHASE** 12* 3X*4HLEG-* 12/ 
1 41H ***** INSUFFICIENT THRUST TO FUi-tOW PATH./////) 

(05 RETURN 

30 IF (SING, GT. 1.0) SING-1. 0 

C CHECK CAPABILITY TO MATCH PATH ALTITUDE. 

DELHC-VEL*0EIV*SING/(32.2*(SING*C2-SING)) 

IF ( ABS ( ( DELHC-DELH) /DELHJ.LE.0.2) GO TU 40 
UO NERROR-2 

WRITE (6*6C2) IPHASE* LEG 

602 FORMATS /////33H *****FATAL ERROR IN PATH. PHASE-* 12* 3X»4HLEG«* 12/ 
1 44 H *****PATH ALTITUDE DEVIATION EXCEEDS LIMIT./////) 

RETURN 
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SUBROUTINE PATHS 76/76 QPT*2 


^.>rir^prr^- 1 .: .. " -t ■ ■■,:•. i ■ — — _s_ 


FTN 4,5+410 05 


15 


2C 


25 


3C 


35 


40 


45 


5C 


55 


6C 


65 


7C 


40 


603 


50 

604 

60 


C“- 


C- 

70 


O 

BC 


go 

60S 

100 


m 

UVR QUALITY 


DELT-DELH/ (VEL+SING) 

CDSG*SQRT£le-SING**2) 

0ELX«VEL*0ELT*CD$G/6Q80» 

DELWF«rN*$FC*DELT/3600. 

TIME-TIME+DELT/60. 

XPATH-XPATH+DELX 
WF-WF+DELWF 
W« W-DELWF 

IF (IPRTdPHASEJ.EQ.O) GO TO 50 

WRITE (6,603) LEG, HI, H2» HN> DELH, SMN1, SHNH, SMN, VEL1, VEL2,VEL,0ELV, 

1 RHO, A; TN,SFC,ALPHA,CDa,KP,Q,Cl,C2,C3»C4,C5,CLl,C6,CL2,CL,CD> 

2 SING, CDSG, DE LHC, DELT,DELX, DEL WF, TIME, XPATH, WF,W 

FORMAT ( /1X,5HLEG -I2/5X,4HH1 »E14.7,5X,4HH2 -E14.7, 5X, 4HHN -E14.7, 
I 3X,6HDELH -E14.7,3X,6HSHN1 »E14. 7/3X, 6HSHN2 -E14.7, 4X, 5HSMH « 

Z E14.7,3X»bHVEll -EI4. 7, 3X> 6HVEL2 -E14 .7,4X,5HVEL -E14.7/ 

3 3X, 6HD ELV -E14 . 7, 4X, 5HRH0 -E14 .7, 6X, 3 HA »El4 . 7, 5X> 4HTN -E14.7, 

4 4X,5HSFC »E14.7/2X,7HALPHA -E14. 7, 4X, 5HCD0 » E14 .7, 5X, 4HK P -E14.7, 

5 6X,3HQ -E14.7,5X,4HC1 - E14 . 7/5X, 4HC2 -E14.7, 5X,4HC3 -E14.7, 

6 5X,4HC4 -E14.7,5X,4HC5 -E14. 7, 4X, 5HC LI »E14.7/5X,4HC6 -E14.7, 

7 4X, 5HCLH -El 4 . 7, 5X, 4HCL -E 14. 7, 5X, 4HC D «E14.7,3X,6HSING -E14.7/ 

0 3X, 6HC OSG «E14.7,2X,7HDELHC -E14. 7, 3X ,6HDELT -E14.7, 3X,6HDELX « 

9 EI4.7, 2X, 7HDELWF -E 14. 7/ 3X, 6HTI ME »E14.7,2X,7HXPATH -E14. 7, 

1 5X,4HHF - t 1 4. 7,6X, 3HW -E14.7/) 

IF (W.GT.WPl) GO TO 60 

IF (KERRUR.EQ.2) WRITE (6,604) IPHASE, LEG 


FORMAT (//25H W.LE.WPL IN PATH. PHASE-, I2,7H LEG-I2/J 

RETURN 

CONTINUE 

HN-H2 

SMN-SMN2 

VEL-VEL2 

BLTIHE-BL T IKfc + TIME/60 . 

BLRANG-BLRANG+XPATH 
— TEST FOR WEAPONS DROP 

IF ( I U (IPHASE ) . EQ.Ci) GO TO 70 

W-W-WMISS 

IWS-0 

— TEST FOR BOMBS DROP 

IF (IBdPHASU .EO.O) GO TO BO 

W-W-WBOMBS 

IBS-0 

— TEST FOR AMMO DROP 

If (lAdPHASE ) .EQ.C) GO TO 90 

W-W-WAMMUN 

IAS-0 

IF (IPRTdPHASE) .EO.O) GO TO 100 
WRITE (6,6C5 ) W,HN,SMN,VEL 

FORMAT! / 6X, 3HW - E14. 7, 5X, 4HHN -E14.7/4X, 5HSMN -E14.7, 

1 4X, 5HVEL -E14.7/) 

SMNT( IPHASE J-SMN 
HNT ( IPHASE ) -HN 
CLT(1PHASE)*CL 
ALPHA! ( IPHASE ) - ALPHA 
WFT( IPHASE)*WF 
TIMET (IPHASE ) -TIME 
VELT(IP:USfc)«VEL 
SFCT (IPHASE ) -S FC 
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TNT ( IPHASE ) *T N 
CDT ( IPHASE ) “CD 
PLOT £lPHASE)«kLD 
bEELTC IPHASE ) *W 
kART£IPHASE)“WTOT 
QT ( IPHASE ) *U 
SFCUT(IPHASE)“SFCU 
TNUTt IPHASE) “THU 
CDINST£IPrtASE)«CDINSP 
PRTOTT £ IPHASE )*■ PRTOT 
XT£IPHASE)«5tPATH 
RhTURN 
END 



it?! 
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SUBROUTINE TRAJUO 


75/76 0PT*Z 


FTN A.5+A10 



SUBROUTINE TRAJOO 

C OUTPUT OF RESULTS FROM TRAJECTORY ROUTINE 
REAL NZST»N2IT*MZ1G*MENDUR*LUIT 

COMMON’ /TRAJCM/ ALPH A* AR W,C D*COL * CDU, CL* DESLF * DR AG* EN *HN* KP* 

1 RLD,SMN*JMAX*RANGE*SFC*SwlNG*THRUST,TW*ULTLF*a*WAMMUN*W£TANK* 

2 WFUEL,WGraWT,WMISS*WTQT*WPLWT*MSTARTE12)»H$TART£12>*TIM£12)* 

3 P5iGT£12)*TDOTST(12)*NZST(12)*PSIF(12)*TDGTIT(12)*NZIT£12J* 

A w8QMBS*WFEXT*SFCU*THRSTU*CDINSP*PRTQT*WKFUEL*CRMACH*FLTQ*X£12)* 

5 FLLAND,TENDUR* BLR ANG* BL TIME ,TI MTGl * TI M102* WFTO* PNAME1 ( 12 ) * 

6 PNAME2(12),SMMT!12)*HNr(12)*WFT(12)*TlMET(12)*XT£12)*CLT£12)* 

7 CDT(12),ALPHAT(12)*CLIT(12)*CDIT(12)*ALPHIT£12J*WFUSE0,WLAND* 

8 tlT£12)»RL0T£12)*SFCT(12)»TNT(12J*CET£12)*0Y(ld)* 

9 IA0*IPS*ITS*IWS»IBS*NPHAS£*IDY(1A) 

COMMON /TRAJEX/ ALP1GT£ 12 ) * CO INS TU2 ) * CDlGTt 12 ) * CL1GTC12 ) * DECEL* 

1 UELHP,D£LMP*DXCRUS,FLFAC*FRFURE*HEND(12),HGNDUR,HMAXP*HMINP* 

2 hPATHT (50)>M£ND(12)*MENDUR*NT£12)*PRT0TT(12)*RADIT(12J* 

3 RADST(12) *SF*SFCUT£12)*SMMAXP*SMMINP* 

A TNUT{12),VlND{12)*tfELT(12)*tfPATHT£50)*WART(12),WCQMBP*WE£LT(12)» 

6 kFRES*wFTOT* WFTQ1* WFT02* WF TR AP * WKL ANO * XOESC * XGRL AN* FWGMAX*TOL* 

6 VFIXED*WGCALC* WPL* WGTa»0UMY(82 )* IA£ 12J *IAS* IB( 12)*I BREG* IENDUR* 

7 IHISS*IP£12)»IPHASE*IPL0T*IPRlNT*IPRTtl2)*IPSIZE*IPST01,lPST02* 

8 lW(12)*IX(12)*KERROR*MILCaM*NCRUSE*NL£GCL*NLEGCR*NLEGLO,NMIS5, 

9 NPATH* I QUMY £ 6 ) 

DATA CUMB/AHCOMB/ 

DATA LOIT/AHLQIT/ 

WRITE ( 6*6CC J IMISS*WPL 

6CC FORMAT f 1H1* 2AX* 17HTRAJ EC TORY OUTPUT// 

1 2AX*7HMISSIQN*I2*12H (PAYLOAD »*F0.D*AH LB)// 

2 67H PHASE M H CL ALPHA WFUEL TIME 

3 VEL/9X,57HSFC£ I) THRUST(I) CD L/D W WA 

A G/9X* 57HS PC (U ) THRUST£U) CDINST PR 

5X) 

DO 10 IPHASE®1*NPHASE 

WRITE £6*601) PNAHE1 ( IPHASE >* PNAHE2 (IPHASE >* SMNT< IPHASE )* 

1 HNT £ I PHASE 1 * CLT ( IPHASE ) *ALPHAT( I PHASE )*WFT£ I PHASE)* TIMET (IPHASE)* 

2 VELT (I PHASE), SFCTt IPHASE )*TNT£ IPHASE ) *CDTt I PHASE) *RLDT( I PHASE)* 

3 hEELT ( I PHASE )* WARTt I PHASE >*QTE I PHASE) *SFCUT ( IPHASE )*TNUT( IPHASE)* 
A COINS T £ IPHASE ) *PRTOTT( I PH ASE)*XT£I PHASE) 

601 FURMAT(/1X*2AA*FA«2*F11.0*F10.A*F8.2*F10.1*FS.2*F8.0/ 

1 Fl3,2*Fllo0*F10,A*F8.2*F10.1*F8.2* F6.0/ 

2 F13.2»F11.0*F10.A*lBX*F8.2*F8.0i 
1C CONTINUE 

WRITE £6*602) WFT01* WFTOT* WFT02* WFEXT * W FUSED* WFUEL* WFRES* WFTRAP 

602 FORHAmHl* 26 x* 12 HFUEL' SUMHARY// 1 AH TAKEOFF FUEL/ 

1 9X,7HWFT01 *s>F8«0,5X»12HT0TAL FUEL **F8,0/ 

2 9X* 7HWFTQ2 «* F8«0*5X,i2H ' EXTERNAL *>*F8.Q/ 

3 16H MISSION FUEL «* F8. 0* 5X, 12H INTERNAL »*F8.0/ 

A 16H RESERVE FUEL «*F0.O/16H TRAPPED FUEL «,F8.0) 

J*0 

NZ1G-1.0 

TD0T1G*0«Q 

RAD1G-0.0 

PSS-0.0 

DO 30 IPHASE«1,NPHASE 
IF (PNAMEKIPHASE) .NE.C0M8) GO TO 30 
IF (J o NE,0)"G0 TO 20 
WRITE ( 6* 603 ) 
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SUBROUTINE TRAJDO 76/76 0PT*2 FTN A.5+AJ.8 


6C3 FORMAT { ///19X*28HADDITI0NAL COMBAT PARAMETERS/ /13X* 

1 58HC0ND IT IGNS PS NZ TOOT RADIUS ALPHA CL CD) 

J-l 

20 WRITE (6*60A) SMNTC IPHASE ) * PS1GT ( IPHAS E ) * NZ1G* TDOTIG^R AD1G* 

1 ALP1GT! IPHASE )*CL1GT( IPHASE ) pC 01GT( IPHASfc ) j HNT { IPHASE 3* PS S» 

2 NZST! I PHASE)*TDOTST( IPHASE >#R ADS TC IPHASE i , ALPHAT { IPHASE) * 

3 C LT C IPHASE ) * CDTl I PHASE > * PS IT 1 1 PHAS t)*NZIT( I PHASE) *TDOT IT C IPHASE)/ 
A KADIT( IPHASE)* ALPHIT! IPHASE )*CLIT( IPHASE) * CO IT ( IPHASE ) * 

5 CET( IPHASE) 

6CA F0RMATI/3H M»*F5.Z*15H 1 G FLIGHT* F8. 1* F 6. 2* F7.2* F8.0* F7. 2* 

1 F7.3»r8.A/3H H«*F6.0*13H SUST AINED* F9. 1* F6. 2* F7. 2* FB. 0* FT. 2* 

2 F7 »3*F8.4/13X*10HMAX. INST»*F8.1*F6.2*F7«2*FBs0*F7«2*F7.3*F8»A/ 

3 13X* 15HC0MBAT ENERGY. »j._E13. 6) 

30 CONTINUE 

WRITE (6*605) BLTIME* BLRANG* FLTO* FL L AND* XGRL AN* WLAND* WGTQ* W* 

1 MENDUR»HENDUR*TENDUR 

605 FCRMATt //12X* 23H BLOCK TIME «*F7.3*7H HOURS/ 

1 35H BLOCK RANGE »*F?.1*7H N. H./ 

2 35H TAKEOFF FIELD LENGTH! TOTAL RUN) «*F7*0*6H FEET/ 

3 35H LANDING FIELD LENGTH! TOTAL RUN) **F7.0#6H FEET/ 

A 35H LANDING FIELD LENGTHCGROUND RUN) =»F7.0*6H FEET/ 

5 35H WEIGHT FOR LANDING CALCULATION **FB.I)*7H POUNDS/ 

6 35H TAKEOFF WEIGHT »*FB,0*7H POUNDS/ 

7 35H LANDING WEIGHT **F8oQ*7H POUNDS/ 

8 35H ENDURANCE MACH NO. “*F5«1/ 

9 35H ENDURANCE ALTITUDE »*F7.Q*6H FEET/ 

1 35H ENDURANCE TIME «*F7o3*7H HOURS) 

DO AO IPHASE »1*NPHASE 
IF IPNANEl(IPHASE).NE.LQIT) GO TO AO 
WRITE ( 6* 6G6) X(IPHASE) 

606 FORMAT (20X*15HL0ITER RADIUS «*F6.2*BH N. M « ) 

AC CONTINUE 

RETURN 





SUBROUTINE TRPLOT 


76/76 OPT® 2 


SUBROUTINE TRPLOTI NERROR, IGEO,KGPRNT > 

REAL LIFT,NZS,NZI 

COMMON /TRAJCM/ ALPHA, ARW, CD, CDL, CDO, CL, DESLF, DR AG, EN, KNjr K P, 

1 RLD,SMN, JHAX, RANGE, S FC,S WING,THRUST,TW, ULTLF » W, WAMKUN, WETANK, 

2 WFUEL,WGTDWT,WHISS,WTQT,WPLWT, MST4RT(121,H5TART(12)»TIH(12), 

3 PSlGT(l2},TDOTSr(12>,NZSTC12J,P5lT{12)»TDOTIT(12),NZIT(lZ), 

A W80MB5 , WFEXT, SFCU, THRSTU, CDINSP, PRTQT,WKFUEL,CRMACH,FLT0,X(12), 

5 FLLAND,TENDUR»BLRANG,BLTIME,TIMT01,TIMTQ2,WFT0, PNAHEH12), 

6 PNAME2 (12 ) , SI1NT( 12 ) ,HNT(1Z) , WFT(12),TIHETI12),XTfl2),CLT(12J, 

7 CDTU2), ALPHATtl2},CLIT(I2),CDIT(12),ALPHIT{I2), WFUSED, WLAND, 

8 QT(123,RLDTC12),SFCTS12),TNT(12),CtT( 123,DY(16), 

9 IA0,IPS,ITS,IWS,IBS,MPHASE,IDY(14) 

COMMON / TRAJEX/ At, PIG T{ 12 1 , COIN ST (12),C01GT(12),CLlGT(12i,DECEL, 

1 DELHP, DELMP, DXCRUS, FLFAC,FRFURE, MEND ( 12 J ,HENDUR, HMAXP , HHINP, 

2 HPATHT(503,MENDC12),MENDUR,NTC121,PRT0TT<12J,RADIT(12), 

3 RAD5T(12>#SF,SFCUT(12),SM(UXP,SHMINP, 

A TNUni2),VINDUZ),VElTClZI,VPAYHT{50),KART(12) , WCOMBP, WEELTUZ), 

5 WFRES, WFTOT* WFTOX, WFT02, WFTRAP, WKL AND, XOESC, XGRLAN, F WGMAX,TOL, 

6 WFIXE0,rfGCALC»WPLiWGIDpDUMr(82»,IAC12}#IAS, IB l 12), IBREG, I ENDUR, 

7 IMISS,IP{12),IPHASE>IPL0T»IPRINT»IPRT(12)»IPSIZE»IPSTD1> IPST02, 

8 IW(12),IX(12),KERR0R,MILC0H,NCRUSE,NLEGCL#NLEGCR,NLEGL0,NMISS, 

9 NPATH, IUUMY £6 ) 

DIMENSION ANS(A) 

DATA CQMB/4HC0MB/ 


— IPLOT * 0, NO PRINT, NO PLOT _ 

— IPLOT = 1, PRINT, NO PLOT 

— IPLOT * 2, NO PRINT, PLOT 

—IPLOT « 3, PRINT, PLOT 

— WCOMBP .LE* 0.0, W " W FROM LAST COMBAT PHASE 
— WCOMBP BETWEEN 0.0 AND 1.0, W ® WC0M6P*WGT0 

— WCOMBP « GT o 1.0, W « WCOMBP 

“•“CALCULATE COMBAT PARAMETERS. 

4=0 

DO 10 IPHASE 3 1,NPHASE 

j®iPHjsE £ltIPHASE,aNe ‘ C ° M) G ° T ° 13 ORIGINAL PAGE IS 

continue OF POOR QUALITY 

IF ( J.EQ.O) I PS“2 
IF (J.GT.3) 1PS-IPCJJ 
JM1® J-l 

IF (JM1.E0.0) JM1»1 

IF ( (WCQMBP.LE.O.O) .AND. i J.EQ.Q) ) W«.7*WGTQ 
IF HWCOMBP.LE.O.OJ.AND.U.GT.On W“WE ELT ( JM1 ) 

IF t (WCOMBP. GT. 6.0) .AND. (WCOMBP. LE, 1.0) J W«WC OMBP*WGTO 
IF (WCOMBP. GT. 1.0) W«WCOHBP 
NH®( NHAXP-HH1NP ) JDELWP+1. 0 
NP* (SHMAXP-SMMINP ) /OELHP+1. 0 

IF ((IPLOT. EQ.IKOR. (IPLOT. EQ. 3)) WRITE (6,600) HMINP, HMAXP, 

1 DELHP, SKMINP, SMMAXP, DEL HP, WCOMBP, IPLOT, QM AX, W 
FORMAT ( 1H1,9H HMINP »F?.0/10H HKAXP *F7»0/10H DELHP -F7.0/ 

1 10H SMMINP *F6. 2/10H SfcHAXP -F6.2/10H DELMP ®F6.2/ 

2 1GH WCOMBP *F9.2/I0H IPLOT »I2/10H QHAX «F7.0/ 

3 10H W a F8 .1//05H H M PS1G NZS TOUT: 

9 NZI PSI TDOTI Q/> 

IF (IPLOT. GE. 2) WRITE (62,601) NH,NH, DELHP, DEIMP, QMAX 
FCiRMAT(21B,3F0.2) 

DO AO IH“I,NH 


ORIGINAL PAGE IS 
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W«.7*WGTQ 
W“WE ELT ( JM1 ) 
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SUBROUTINE TRPLQT 


7b/7fc QPT*2 


FTN 4.5+410 


HMFLOATt IH-1)*0ELHP 
CALL AT 62 (HNj ANS J 
RHO-ANSQ) 

A*AN5 f 4 ) 

OU 4C 

SKN-FLQATdh-D+DELHP+SHKINP 

V£L«SMN*A 

Q-.5*RHQ*VEL**2 

QS*Q*SklNG 

CL»N/OS 

IA0*B 

CALL TRAU01I2) 

CALL AER0(2,N£RRQRj IGEOn KGPRNT) 
IF (NERRQR.GE.2) RETURN 
CALL TRAJOl(l) 

DR AG»CD*QS 
CALL TRAJOI 1 2 ) 

CALL P RQP ( 2* MERRQRjIGED* KGPRNT) 
IF (NERR0R.GE.2) RETURN 
CALL TRAJ01C1) 

TN»THRUST*cN 

PS1G*(TN-DRAG)*VEL/W 

N2S-1.001 

IF (PS1G.LE.0.) GO TO 20 
I AQ“4 

CALL TRAJ0H2) 

CALL A£R0(2>NERRQR>IGE0» KGPRNT) 
IF { NERROR «GE .2) RETURN 
CALL TRAJ31U) 

LIFT*CL*tfS 

NZS“ (LIFT+TN+SIN (ALPHA/57. 3) ) /W 
IF (NZS.LE.DESLF) GO TO 20 
IAO-6 

CALL TRAJ01C2) 

CALL AEROC 2>NERRQR> IGEO»KGPRNT) 
IF (NERROK o GE • 2} RETURN 
CALL TRAJOl(l) 

NZS*DESLF 

20 IF (NZS.LE.I.O) NZS-1.001 
PHI«ARCQS(1.0/NZS) 
RADIUS»VEL**2/<32.2*TAN{PHI]) 
TD0TS*57.3*V£I /RADIUS 
IA0=5 

CALL TRAJ0K2) 

CALL AERO ( 2* NERROR; I GEO# KGPRNT) 
IF (NERR0R.CE.2) RETURN 
CALL TRAJOI(l) 

LIFT-CL+QS 

NZI- (LIFT+TMSIN I ALPHA/ 57. 3) )/U 
IF (NZI.LE.DtSLF) GO TO 3C 
I A0*6 

CALL TRAJ0K2) 

CALL AER0(2#NERR OR >IGEOi> KGPRNT) 
IF (NERR0R.GE.2! RETURN 
CALL TRAJOltl) 

NZI*DESLF 

IF (NZI.LE.1.0) NZI-I.001 
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SUBROUTINE TRPLOT 


76/76 QPT-2 


FTN A. 5+410 


05 


15 


2C 


25 


30 


35 


AC 


AS 


5C 


55 


602 

AO 


DRAG«vO*OS 

PSI* ( TN*CQS (ALPHA/57.3)-DRAG)*VEL/W 
PHI-ARCOSd.O/NZI ) 

RADIUS -VEL**2/< 32, 2*TAN( PHI) ) 

TD0TI-57.3+VEL/RADIUS 

IF UIPL0T.EQ.D.0R.CIPL0T.EQ.3)) WRITE (6,602) HN,SMN, 

1 P$1G,NZS,TD0TS,NZI,PSI,TD0TI,Q 
IF (IPL0T.GE.2) WRITE (42,602) HN,SHN, P51G,NZS,TD0TS,NZI, PSI, 

1 TUOTI, Q 

FORMAT (F7.0,F6»2, 1P7E11.3) 

CONTINUE 

C LANDING FIELD LENGTHS AT ALTITUDES OF 0, 4000, 8000, AND 12000 FT. 

IF UlPLOr.EQ.l>.0{MIPLOT.EQ.3)) WRITE (6,603) 

603 FORMAT <1H1,8X,33HLANDING FIELD LENGTH CALCULATIONS// 

1 tcH ALT „ WMIN FLLAN(WHIN) WGTO FLLANt W 

26T0) // ) 

WhlN-WGTO-WFTOT-WPL 

IAD-12 

DO 50 1-1,4 . . . .. 

HN»(I-1)*4000 

CALL AT 62( HN, ANS ) __ 

RHO- ANS ( 1 ) ORIGINAL PAGP TQ 

CALL TRAJQK2) OP POOR ^7^ 

CALL AERO(2,NERROR,IGEO, KGPRNT) QUALITY 

IF (NERR0R.GE.2) RETURN 

CALL TRAJOl(l) 

VSTALL-SORl ( 2 . *WMIN/ ( SWING+QL*RHQ ) ) 

VSCREN-1.3 + VS TALL 
VTDOWN-1, 15*VSTALL 

XAIR«({VSCREN**2-VTDQWN**2)/64.4+5D.>*RLD 
XGRLAN-VT0QWN**2/ ( 64* 4<*DECEL ) 

FLLAN1 » ( XA1 R+XGRLAN) /FLFAC 
VSTALL«S0RT(2.*WGTD/(SWING*CL*RH0) ) 

VSCREN«1.3*VSTALL 

VT00WN«1.15*VSTALL 

X AIR- ((VSCREN++2-VTD0WN++2) /6A. A + 50. )+RLD 
XGRLAN-VT00WN++2/ (6A.A4DECEL) 

FLLAN2- (XAIR+XGRLAN) /FLFAC 

IF ( (IPLOT.EQ.l) .QR.UPL0T.EQ.3)) WRITE (6,60A) HN,WMIN,FLLAN1, 

1 WGT0,FLLAN2 
6C A F0RMAT(F7.C,AE15e6) 

IF (IPL0T.GE.2) WRITE (42,604) HN,WMIN,FLLAN1,WGT0, FLLAN2 
50 CONTINUE 
RETURN 
ENO 
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SUBROUTINE GOLDEN 76/76 OPT-2 FTN 4.5+410 < 


1 SUBROUTINE GOL DEN t XMIN, XHAX, I XY , NMA X , X X, Y, NE RROR , IGEO, KGP RNT ) 

C THIS ROUTINE USES THE GOLDEN SECTION METHOD TO FIND 

C THE X CORRESPONDING TO THE MAXIMUM Y IN THE INTERVAL 

C BETWEEN XMIN AND XMAX. 

5 C IF IXY - 1# THEN X * ALT. AND Y - BREG. FACTOR 

C IF IXY - 2, THEN X » ALT. AND Y - ENDUR. FACTOR 

C IF IXY ■ 3, THEN X « MACH AND Y • ENDUR. FACTOR 

COMMON /TRAJCM/ AL PHA, AR W, CO# CDL, CDQ, CL, DE S L F * DR AG, E N, HN, KP, 

1 RLD,SMN,QMAX, RANGE ,SFC, SWING. THRUST# TW,ULTLF>W,WAHMUN,WETANK, 

10 2 WFUEL,WGT0WT,WHISS,WT0T,WPLWT,HSTART(12),HSTART[12) >TIM(12), 

3 PS1GT(12),TD0TST(12),NZST(12),P5IT(I2), TDQTITf 12),NZIT(12), 

A WBOHBS.WFEXT, $FCU, THRSTUoCDINS P, PRTOT, WRPUfc L, CRHACH , F LTO , X { 12 ) , 

5 FLLAND,TENDUR,BLRANG,BLTIME, TI HT01 , T I HT02, WF TO, PN AME1 l 12 ) , 

6 PNAME2(12)»SMNT(12J,HNr(12) ,WFT{12) , TIMET ( 1 2 ) , XT (12 ) ,CLT ( 12 >, 

15 7 CDT(12),ALPHAT(12),CLIT( 12) ,CUIT£12),ALPHIT( 12) ,WFUSED,WLAND, 

8 QT(12),RL0T(12),SFCT(12),TNT(12),CET(12),0Y(18)» 

9 IA0,IPS,ITS,IWS,IBS,NPHASE, IDY 1 14 ) 

COMMON /TRAJEX/ AL P1GT ( 12 ) , CD I NS T 1 12 >, CD1GT f 1 2 ) » CL1GT ( 12 1 , DECEL, 

1 DtLHP, DfcLMP,DXCRU5,FLPAC, FRFUR E , HEND ( 12 ) , HE NDUR i HMAXP, HM INP , 

2C 2 HPATHT { 50 J , MEND ( 12 ) , MENDUR, NT ( 12 ) , PRTOTT (12)»RA0IT(12), 

3 RADSTt 12) , SF,SFCUT( 12) , SMNAXP, SMMINP, 

4 TNUTU2). VlNDt l2)>VELT(12)f VPATHTt 50 ) , WAR T { 12 ) , WCOMB P, WE ELT (1 2 1 , 

5 WFRES,WFTQT,WFT01,WFT02,WFTRAP,WK LAND, XDESC/XGR LAN, FWGMAX.TOL, 

6 WFIXED,wGCALC,WPL>WGT0,DUHY(82 ) > I A ( 12 ) , I AS, I B ( 12 ) , I BREG, I ENDUR# 

25 7 I MISS, IP ( 12 KIPHASE, IP LOT, IPRINT.IPRT (1Z),IPSIZE, IP ST 01, IPST02, 

B IW(12),IX(12),ICERROR,MILCOM,NCRUSE#NLEGCL,NLEGCR,NIEGLO,NMISS, 

9 NPATH, IDUMYS6) 

DIMENSION ANS !4J 
IPSA VE*IP$ 

30 KCALC-2 

IAQ« 9 
I PS* A 

IF (fIPRINT.GE.2KAND.tIXY.EQ.in WRITE (6,600) 

IF UIPRINT.GE.2).AND.(IXY.EQ.2) ) WRITE (6.601) 

35 IF ((1PRINT.GE.2). AND .(IXY.EQ.3)) WRITE (6,602) 

600 FORMAT ( /22H SEARCH FOR BREG. ALT./23H ALT. BREG. FACT., 

1 5X,4HHACH,8X # 2HCL>UX,2HCD, 7X, 6HTHRUS T , 3X, 3HS FC ) 

601 FORMAT ( /23H SEARCH FOR ENDUR. ALT./23H ALT. ENOUR. FACT., 

1 5X»4BMACH,8X , 2HCL, 8X, 2HCD, 7X, 6HTHR UST, 3X, 3HS FC ) 

AC 602 FORMAT ( /23H SEARCH FOR ENDUR. HACH/23H MACH ENOUR. FACT., 

1 5X, AhMACH,0X,2HCL,8X,2HCD,7X,6HTHRUST,3X, 3HSFC ) 

X 1*XMIN 
X2*XHAX 
XX«X1 

AS N*0 

XBbST-XX 

YBEST—i. 

IF (IXY.LE.2) HN-XX 
IF (IXY.EQ.3) SMN*XX 
50 IRETRN«1 

C START OF FUNCTION EVALUATION 

C 

10 N-N+l 

55 CALL AT62 (HN,ANS ) 

RHO-ANS (1) 

A«ANS ( A ) 
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SUBROUTINE GOLDEN 76/76 0PT*2 


FIN 4.5+410 01 


60 


65 


70 


75 


ac 


85 


90 


95 


oc 


05 


.1C 


603 


C- 

c 

c- 

20 


3C 


4 0 
50 


6C 


7C 

8C 


VEL-SMN+A 

CL*2#*W/(RH0*SWING*VEL**2) 

CALL TRAJ0K2) 

CALL AERO(KCALC?NERROR# I GE 0> KGP RNT ) 

IF tNERRDR.GE.2) RETURN 
CALL TRAJ01C1J 
DRAG=«5*RHD*SWING*CD*VEL**2 
CALL TRAJ0K2) 

CALL PRJP{KCALC*NERRORj. IGEO jKGPRNF) 

IF t NERROR . GE .2 ) RETURN 
CALL TRAJOKU 

IF (IXY.EO.l) Y-3600.+VEL+RL0/SFC 
IF CIXY.GE.Z) r-RLD/SFC 
IF (NERROR.EQ.il Y—1.0 

IF ( 1PRINT.GE.2) WRITE {6^603) XX# Y# SHN# CL # C Dj. THRUST# SFC 
FCRNAK lX#G9e3#G13.4#G14.3#2Gl2.4#G13.6#G12. 3) 

IF ( Y.GT. YBEST) XBEST-XX 
IF ( Y. GT * YBES T) YBEST* Y 


END OF FUNCTION EVALUATION 


GP TO (20#30#40#6Q)# IRETRN 

Y1*Y 

XX*X2 

IF ( IXY • LE « 2 ) HN*XX 
IF (IXY.EQ.3) SMN-XX 
IRETRN* 2 
GO TO 10 
Y2 * Y 

TAU*.5*(SQRT(5. ) -1 * ) 
X3*X2-TAU*(X2~X1) 

X X *X3 

IF (IXY.LE.2) HN*XX 

IF (IXY.EQ.3) SMN* XX 

IRETRN-3 

GO TO 10 

Y3-Y 

X12*.5MX1 + X2) 

DELX*ABS(X12-X3) 

X4-X1Z+DELX 

IF U3.GT.X12) X4“ X12-DELX 
XX*X4 

IF (IXYdLE.2) HN-XX 

IF (IXY.EQ.3) SNN»XX 

IRETRN-4 

GO TO 10 

Y4-Y 

IF (Y3.GE.Y4) GO TO 90 

IF (X4.LT.X3) GO TO 70 

X1*X3 

Y1*Y3 

GO TO ao 

X2*X3 
Y2-Y3 
X3-X4 
Y3*Y4 
GO TO 110 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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SUBROUTINE GOLDEN 76/76 UPT”2 


FTN A.S+'tlO 


IF (X3.LT.XA) GO TO 100 

Xl-XA 

Y1*YA 

GO TO 110 

X2-X4 

Y2*YA 

IF 1N.LE.NMAX) GO TO 50 

XX-XBEST 

Y*YBEST 

IF (Y.EQ.-l.O) NERRQR-2 

IF (NERROR.EQ.23 WRITE (6j60A) XX»Y>IXY 
FORMAT ( 35H ***FATAL ERROR IN GOLDEN ROUTINE./ 

1 5X,3HX =>E12.5>2X,3HY E 12 . 5# ZX, 5HI X Y *>I5) 

IF niPRINT.GE.2J. OR. (KERROR.EQ. 21) WRITE (6*605) XX*Y 
FORMAT C / 35H MACH OR ALTITUDE SEARCH CONCLUOED./ 

1 5X* BHX «,E12.5*2X*3HY -*E12.5! 

IPS«IPSAtfE 
RETURN 
END 
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